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ABSTRACT 


This  report  describes  an  architecture  for  a  highly  parallel  system  of  compu¬ 
tational  processors  specialized  for  real-time  adaptive  antenna  nulling  computations 
with  many  degrees  of  freedom,  called  MUSE,  and  a  specific  realization  of  MUSE 
for  64  degrees  of  freedom.  Each  processor  uses  the  CORDIC  algorithm  and  has  been 
designed  as  a  single  integrated  circuit.  Ninety-six  such  processors  working  together 
can  update  the  64-element  nulling  weights  based  on  300  new  observations  in  only 
6.7  ms.  This  is  equivalent  to  2.88  Giga-ops  for  a  conventional  processor.  The  com¬ 
putations  are  accurate  enough  to  support  50  dB  of  signal-to-noise  improvement  in 
a  sidelobe  canceller.  The  connectivity  between  processors  is  quite  simple  and  per¬ 
mits  MUSE  to  be  realized  on  a  single  large  wafer,  using  restructurable  VLSI.  The 
complete  design  of  such  a  wafer  is  described. 
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1.  INTRODUCTION 


When  an  array  of  N  antenna  elements  is  subject  to  undesired  interference,  such  a s  co-channel 
interference  or  jamming  plus  the  thermal  noise  in  each  of  the  N  receivers,  the  interference  power 
can  be  reduced,  relative  to  the  power  in  some  desired  signal,  by  forming  a  suitable  weighted  sum  of 
the  waveforms  observed  on  all  the  antenna  elements  as  the  system  output.  This  is  called  nulling. 
In  many  instances  the  choice  of  suitable  weights  must  be  made  adaptively  with  time.  The  optimum 
choice  of  weights,  in  the  sense  that  the  signal-to-interference  ratio  observed  in  the  system  output 
is  maximized,  is  the  solution  to  a  well  studied  least-squares  problem.  The  number  of  arithmetic 
steps  required  to  solve  this  least-squares  problem,  almost  independent  of  the  algorithm  chosen, 
is  proportional  to  the  cube  of  the  number  of  antenna  elements.  Furthermore,  in  the  presence  of 
interference  that  may  change  with  time,  it  is  necessary  to  perform  the  adaptive  weight  determination 
repetitively  and  in  real  time. 

When  the  antenna  is  on  board  a  moving  spacecraft,  the  combination  of  a  real-time  requirement 
driven  by  the  satellite  motion  and  the  cubic  dependence  of  the  computational  cost  of  adaptive  weight 
determination  gives  a  limit  to  the  number  of  antenna  elements  that  may  practically  be  nulled.  A 
previous  study  of  this  computational  cost  [1]  set  this  limit  at  about  N  =  26,  based  on  an  assumed 
conventional  digital  signal  processing  architecture.  Such  a  limitation  is  in  no  sense  absolute,  since 
it  depends  on  what  resources  we  are  willing  to  allocate  to  a  nulling  processor.  However,  because  the 
difficulty  of  the  nulling  computation  grows  as  iV3,  simply  using  more  resources  is  not  an  efficient 
way  to  handle  a  large  number  of  antennas. 

In  this  technical  report,  a  specialized  adaptive  nulling  processor  called  MUSE,1  which  is  capable 
of  determining  the  weights  for  N  =  64  antenna  elements,  is  described.  Due  to  its  novel  architecture, 
it  may  be  realized  compactly  using  restructurable  wafer  scale  integration;  it  would  then  be  the  size 
of  the  square  inscribed  in  a  5-in-diam.  circle.  MUSE  is  not  only  substantially  smaller  and  lighter 
than  a  conventional  processor,  but  it  also  will  use  substantially  less  power.  Although  the  MUSE 
processor  is  specialized  for  N  =  64  antennas,  the  design  concept  of  MUSE  can  be  applied  to  design 
a  similar  processor  for  an  antenna  array  with  either  more  or  fewer  elements. 

In  the  next  section  there  is  a  brief  explanation  of  the  mathematics  of  adaptive  nulling  and 
the  method  of  voltage  domain  computation  of  a  Cholesky  factor  using  Givens  transformations.  In 
section  3,  there  is  an  explanation  of  CORDIC  rotations  and  a  demonstration  of  how  they  are  suitable 
for  realizing  the  Givens  transformations  needed  in  updating  a  Cholesky  factor.  In  section  4,  the 
idea  of  a  “naive”  systolic  array  of  computing  elements,  each  of  which  is  composed  of  three  CORDIC 
rotation  cells  operating  together,  is  developed.  So  connected,  an  array  of  computing  elements 
is  capable  of  sharing  the  work  of  updating  a  Cholesky  factor,  which  is  the  largest  part  of  the 
computational  task.  However,  it  is  relatively  inefficient.  In  section  5,  we  show  how  to  modify  the 
systolic  array  to  make  it  100  percent  efficient  for  the  Cholesky  update  task,  using  a  technique  which 
is  called  latency-controlled  interleaving.  In  section  6,  the  parameter  choices  for  the  CORDIC  cell, 
influenced  by  the  specific  system  goals,  are  explained.  A  simulation  used  to  motivate  the  choice  of 

1  MUSE  is  an  acronym  which  stands  for  Matrix  Update  Systolic  Experiment. 
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some  of  these  parameters  is  described.  In  section  7,  we  describe  some  of  the  technology  issues  of 
restructurable  VLSI  and  show  how  the  systolic  array  is  suited  to  implementation  in  this  technology. 

In  section  8,  the  problem  of  solving  the  linear  equations  for  the  optimum  weights,  once  a 
Cholesky  factor  is  found,  is  discussed.  CORDIC  cells  can  be  used  for  this  computation  in  a  manner 
that  is  quite  closely  related  to  the  way  in  which  they  are  used  to  update  a  Cholesky  factor.  A 
modification  is  described  of  the  timing  of  the  systolic  array  so  that  it  may  be  used  for  both  the 
Cholesky  update  and  the  solution  of  the  linear  equations. 

In  section  9,  a  detailed  description  is  provided  of  a  cell  that  performs  the  CORDIC  function 
and  all  related  necessary  functions  for  building  up  a  MUSE  processor  on  a  wafer  from  96  identical 
cells  of  that  type.  In  section  10,  some  specifics  are  given  of  the  topology  of  such  a  wafer  and  of 
the  method  by  which  the  cells  are  connected  together  to  form  a  working  system.  In  section  11, 
the  architecture  of  a  test  and  demonstration  system  for  demonstrating  the  performance  of  such  an 
integrated  wafer  is  described. 
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2.  VOLTAGE  DOMAIN  NULLING  WITH  GIVENS  TRANSFORMATIONS 


Figure  2-1  is  a  cartoon  of  a  typical  adaptive  nulling  system.  The  N  antenna  elements  each 
produce  a  waveform  that  is  assumed  to  be  down-converted,  sampled,  and  digitized.  At  the  n-th 
sampling  instant  a  complex  number  Xin  comes  from  the  ith  antenna.  All  the  antennas  are  assumed 
to  be  sampled  at  the  same  instant,  so  we  may  collect  together  the  samples  from  the  same  sampling 
instant  into  a  column  vector  X_(n).  These  vector  samples  are  passed  to  the  right,  where  the  nulling 
takes  place,  and  some  samples  are  also  passed  down  to  where  the  adaptation  takes  place. 


144248-1 


Figure  2-1.  Cartoon  of  typical  adaptive  nulling  system. 

The  adapted  weights  are  also  collected  into  a  column  vector,  which  will  be  labeled  as  W. 
Although  W  will  be  time- varying,  it  will  not  be  recomputed  for  every  sampling  interval.  A  block  is 
defined  as  the  interval  within  which  W  will  be  held  constant.  The  nulled  output,  y(n),  a  scalar,  is 
formed  by  using  weights  W_  on  the  vector  A(n).  We  can  use  the  vector  notation 

y(n)  =  WhX(n), 


where  ()h  means  Hermitian  transpose.  This  notation  implies  that  we  actually  multiply  the  elements 
of  2L(n)  by  the  complex  conjugates  of  the  components  of  IT. 
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In  Figure  2-1  the  adaptation  process  involves  samples  X(n)  and  also  a  lower  triangular  matrix, 
Ljt,  which  is  regularly  updated  and  then  used  to  determine  W.  The  process  of  adaptation  will  now 
be  explained. 

Within  each  block,  the  optimization  of  W_  is  with  respect  to  the  statistics  of  interference 
appropriate  to  that  block.  The  statistics  that  matter  are  the  correlations  of  interference  observed 
from  one  antenna  element  to  another,  arranged  into  a  matrix  R.  A  mathematical  description  of 
the  nulling  system’s  performance  will  refer  to  this  correlation  matrix  of  the  interference,  which  is 
generally  not  known  in  practice  since  it  is  a  statistical  expectation. 

R  =  £  (j(n)Jh{n)) 

where  J(n)  is  the  interference  component  of  2C(n).  R  is  never  perfectly  known  for  two  reasons: 

•  One  sees  the  interference  and  the  desired  signal  at  the  same  time  as  an  insep¬ 
arable  sum. 

•  One  has  neither  the  time  to  gather  enough  samples  nor  the  computational 
resources  to  use  all  the  samples  available  to  estimate  R  perfectly. 

We  are  going  to  make  the  simplifying  assumption  that  although  the  observed  antenna  signals 
X(n )  contain  both  desired  signals  and  interference,  the  desired  signals  are  very  much  weaker  than 
the  interference  and  may  be  neglected  in  forming  the  correlation  matrix.  This  simplifying  assump¬ 
tion  is  reasonable  because  a  radar  designer  always  tries  to  get  by  with  the  least  possible  transmitter 
power  compatible  with  detecting  a  target  after  signal  processing.  Since  signal  processing  enhances 
the  target  return  relative  to  interference,  it  follows  that  the  designer  will  have  chosen  the  transmit¬ 
ter  power  small  enough  that  the  target  return  before  any  signal  processing  will  be  below  the  noise. 
Therefore,  the  first  limitation  is  not  a  problem1  and  only  the  second  is  serious.  A  limited  amount  of 
data  must  be  used  to  obtain  an  estimate  of  R  for  each  block  and  then,  since  R  is  changing,  another 
limited  set  of  data  may  be  used  each  time  R  is  updated,  and  so  on. 

Let  us  therefore  collect  together  all  the  samples  of  X_(n)  that  we  intend  to  use  in  determination 
of  “optimum”  weights  for  a  given  block  and  refer  to  this  collection  of  data  as  X .  There  is  no 
implication  that  all  the  samples  used  come  from  within  the  given  block  -  we  actually  envision  using 
samples  from  within  a  sliding  window.  If  the  number  of  samples  used  is  M  and  each  sample  is  an 
N  element  vector,  it  is  natural  to  arrange  the  M  vectors  into  a  matrix  of  M  columns  and  N  rows. 
The  estimate  of  the  correlation  matrix  based  on  this  matrix  of  raw  data  is 

R  =  ^-rXXh. 

M 

1  This  argument  assumes  that  only  small  targets  are  to  be  detected.  If  a  mix  of  small  and 
large  targets  are  sought,  it  may  not  be  appropriate  to  neglect  the  largest  targets  when  estimating 
interference  statistics. 
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Since  we  want  to  optimize  a  signal-to-noise  (S/N)  ratio,  we  need  to  characterize  the  signal. 
Let  S  be  a  known  constant  vector  related  to  the  desired  signal  as  follows:  if  all  interference  were 
absent  and  only  signal  present,  we  would  expect  to  observe,  in  X[n),  some  scalar  time  function 
times  the  vector  5. 

Then  the  optimum  choice  of  W_  is  known  from  the  literature  [2]  to  satisfy  the  equation 
RW  =  S. 

This  equation  is  of  great  importance  to  us  for  two  reasons.  First,  the  only  information  needed 
from  the  block  X  to  determine  the  optimum  weight  vector  is  that  information  Avhich  could  be  used 
to  estimate  R.  Hence,  if  we  intend  to  estimate  R  using  R  =  j,xx\  we  can  be  comfortable  in 
transforming  X  in  any  way  so  long  as  that  transformation  leaves  XXh  unchanged.  Second,  the 
equation  tells  us  that  we  can  find  IT  from  R  by  solving  linear  equations  and  hence  it  suggests  that 
we  might  accept  as  a  pretty  good  estimate  [3]  of  W  the  solution  to  the  same  equations  using  R, 

RW  =  S. 

How  can  X  be  transformed  while  leaving  XXh  unchanged?  Suppose  we  postmultiply  X  by  a 
unitary  matrix  Q : 

X  =  XQ. 

By  definition  a  unitary  matrix  has  the  property 

QQh  =  I 


so  we  may  insert  QQh  between  X  and  Xh  in  the  equation  for  R 
R  =  XXh  =  XQQhXh  =  XXh. 

That  means  X  can  be  transformed  into  X  without  changing  the  correlation  matrix  that  is  estimated 
from  it.  Furthermore,  this  transformation  may  be  repeated  again  with  another  unitary  matrix,  and 
still  again,  as  often  as  liked,  until  the  final  transformed  version  of  X  has  a  convenient  form. 

The  strategy  will  be  to  make  successive  choices  of  Q  that  cause  selected  elements  in  the  matrix 
XQ  to  become  zeroes,  but  such  that  elements  of  X  that  have  already  been  zeroed  by  previous 
transformations  are  not  changed.  A  series  of  such  transformations  can  confine  the  nonzero  elements 
of  the  N  x  M  matrix  X  to  an  JV  x  JV  lower-triangular  submatrix,  which  is  called  L,  on  the  left, 

e.g., 

XXX 
XXX 
XXX 
XXX 


X 

X 

X 

X 


X 

X 

X 

X 


X  X 
X  X 
X  X 
X  X 


X 

X 

X 

X 


X 

X 

X 

X 


X 

X 

X 

X 


=> 

■  z 
z  z 

z  z  z 

o 

_  l  l  l  l 

=  [L|0] 
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The  matrices  Lk- 1  and  L *.  in  Figure  2-1  are  instances  of  L. 

Having  made  this  sequence  of  transformations  we  can  write 

i  xx'-A-±ll\ 

and  then  the  linear  equations  we  must  solve  for  W_  take  the  simpler  form 
LLhW  =  MS. 

The  known  constant  M  may  be  absorbed  into  the  steering  vector  S_  and  is  hereafter  omitted.  (Note 
that  since  W_  is  the  solution  to  a  set  of  linear  equations,  if  5  is  scaled  by  any  constant,  we  scale 
the  solution  W_  by  the  same  constant.  However,  the  signal- to-interference  ratio  in  WhX(n)  is 
unchanged  by  multiplying  W_  by  a  scale  factor.) 

We  can  solve  LLh\V  =  5  in  two  easy  steps.  First  define  an  intermediate  vector  variable  Y 
defined  by 

Lh\V  =  Y ,  (2.1) 

so  that  Y  is  the  solution  to 

LY  =  S.  (2.2) 

We  can  first  solve  Equation  (2.2)  for  T  and  then,  having  found  T,  we  are  in  position  to  solve 
Equation  (2.1)  for  W_.  Both  these  equations  are  easier  to  solve  than  general  matrix/vector  equations 
because  the  matrices  involved  are  triangular.  If  we  had  to  solve  a  matrix/vector  equation  involving  a 
general  TV  x  TV  square  matrix,  we  would  expect  to  require  a  number  of  arithmetic  steps  proportional 
to  N 3.  However,  when  the  matrix  is  triangular,  the  straightforward  method  of  solution  requires  a 
number  of  arithmetic  steps  proportional  to  only  TV2. 

To  illustrate  this,  consider  a  4  x  4  special  case  of  Equation  (2.2)  and  Equation  (2.1): 


hiVi  = 

Sl 

(2.3) 

hlVl  +  ^22  2/2  = 

S2 

(2.4) 

hlVl  +  ^322/2  +  ^332/3  = 

S3 

(2.5) 

^41 2/1  +  ^422/2  +  ^432/3  +  ^442/4  = 

S4 

(2.6) 

I44W4 

=  2/4 

(2.7) 

I43W4  +  I33W3 

=  2/3 

(2.8) 

I42W4  ■+■  £32^3  ■(■  £22^2 

=  2/2 

(2.9) 

£4^4  +  £3^3  +  I21W2  +  £n»i 

=  2/1 

(2.10) 

From  Equation  (2.3)  we  immediately  get 
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Si 

y\  =  T 
n  1 


and  the  now  known  value  of  y\  maybe  used  in  the  next  three  equations,  giving: 

h 22/2  =  s2  ~  hiVi  =  s2 

^322/2  +  ^332/3  =  S3  -  ^312/1  =  S3 

^4 22/2  +  ^432/3  +  ^442/4  =  S4  —  ^412/1  =  s4 

This  shows  that  we  solve  for  one  unknown  in  an  N  x  N  triangular  system  and  also  reduce  the 
problem  to  an  (N  —  1)  x  (N  —  1)  triangular  system  at  a  cost  of  N  —  1  multiply-and-add  operations 
and  one  division. 

Now  we  immediately  find  2/2  =  S2/I22  and  reduce  to  a  2  x  2  triangular  system  in  an  analogous 
manner.  Taking  this  process  to  the  end,  we  will  need  to  use  N  divisions  and  (jV—  1)+(N— 2)  +  . .  .  +  1 
multiply-and-add  operations.  The  series  sums  to  | N(N  —  1),  which  approaches  \N2  for  large  N. 

After  using  Equations  (2.3)  through  (2.6)  to  find  2/1 , 2/2 1 2/3 1 2/4 1  these  can  be  used  as  the  known 
right-hand  side  in  Equations  (2.7)  through  (2.10).  Thus 


104  = 


y± 

/*  ’ 
‘44 


and  the  now  known  value  of  104  may  be  used  in  Equations  (2.8),  (2.9),  and  (2.10). 


^32w3  +  l22w3 
/31I03  +  I21W3  +  *nu’l 


2/3  -  *43^4  =  2/3 
2/2  “  *42^4  =  2/2 
2/1  -  *41^4  =  2/i 


Now  we  have 


wz  = 


y%_ 

ih' 


and  we  may  continue  in  an  analogous  manner.  The  operation  count  for  solving  the  second  triangular 
set  of  equations  is  identical  to  that  for  the  first  set.2 

The  above  approach  to  adaptive  nulling  may  be  broken  into  several  phases: 


Obtain  typical  interference:  Form  a  rectangular  matrix  X  whose  columns  are  the 
(vector)  samples  of  interference  observed  at  the  antenna  elements. 

2  If  the  antenna  array  is  designed  such  that  the  appropriate  steering  vector  is  5  =  [0, 0, . . . ,  0, 1]*, 
the  first  triangular  set  of  equations  becomes  trivial  andonly  the  second  set  needs  to  be  solved.  This 
kind  of  adaptive  antenna  array  is  known  as  a  sidelobe  canceller. 
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TViangularize:  Use  orthogonal  transformations  to  compute  a  triangular  matrix  L  from 

X. 

Back-substitute:  Solve  the  set  LF  =  S_  for  the  intermediate  variables  F. 

Back-substitute:  Solve  the  set  LhW  =  F  for  the  desired  weight  vector  W_. 

This  may  be  contrasted  with  a  mathematically  equivalent  approach  in  which  we  compute 

R  =  XXh 

followed  by  a  procedure  called  Cholesky  factorization  for  determining  L  from  R  and  then  carry 
out  the  two  back-substitutions.  The  latter  approach,  while  valid,  is  not  numerically  robust.  In 
the  matrix  R,  the  largest  and  smallest  eigenvalues  are  usually  determined,  respectively,  by  jammer 
power  and  thermal  noise  power.  This  ratio,  if  very  large,  determines  the  minimum  word-length 
needed  in  numerical  operations  to  solve  linear  equations  involving  R.  This  word-length  limitation 
is  fundamental  and  does  not  depend  on  the  algorithm  used  to  solve  the  linear  equations.  But  if  L 
is  determined  directly  from  the  raw  data  A",  without  ever  computing  /?,  the  minimum  word-length 
requirement  for  solving  linear  equations  involving  L  or  Lh  is  about  half  that  needed  with  R. 

Any  nulling  algorithm  that  depends  on  explicit  computation  of  R  from  X  is  called  a  power 
domain  algorithm.  By  contrast,  an  algorithm  that  determines  L  by  transforming  the  data  set  X  is 
called  a  voltage  domain  algorithm.  Because  it  is  most  economical  to  design  and  construct  hardware 
with  short  word-lengths,  voltage  domain  algorithms  are  preferred. 

Now  let  us  give  some  attention  to  the  treatment  of  successive  blocks.  For  this  we  must  recognize 
that  both  X  and  L  will  be  different  in  different  blocks.  Let  the  raw  data  and  the  Cholesky  matrix 
in  the  kth  block  be  X *.  and  L*.  But  how  is  X^  related  to  Xk-i? 

If  the  statistics  of  interference  were  expected  to  remain  constant  with  time,  then  we  would 
expect  to  estimate  those  statistics  better  and  better  by  using  more  and  more  raw  data  in  the 
successive  blocks.  \Ve  would  append  the  observations  that  became  available  during  the  kth  block, 
XneW >  onto  the  matrix  of  all  the  observations  available  before,  X^-i'- 

Xk  =  [Xfc_!  I  Xneu)] 

so  that 

Rk  =  R-k  —  1  "h  X new  X. new 

and 


LkLk  =  Lk~iLlk_1  +  XnewX%ew 
LkLhk  =  [Lfc-j  |  Xnew\\Lk. i  |  Xnew]h. 
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This  suggests  an  efficient  way  to  update  L  by  use  of  unitary  transformations.  If  Q  is  a  unitary 
transformation,  then 


LkLhk  =  [!*_!  | 

Xnew  iep*-i  i  Xnew  ]Q)h- 

We  choose  Q  so  as  to  force  zeroes  into  the  positions  occupied  by  Xnew,  thus 


l 

n 

n 

n 

n 

n 

n 

l 

l 

l 

n 

n 

n 

n 

n 

n 

— K. 

l 

l 

l 

l 

l 

n 

n 

n 

n 

n 

n 

l 

l 

l 

l 

l 

l  l 

n 

n 

n 

n 

n 

n 

l 

l 

l  l 

There  are  two  advantages  to  updating  L*  using  Lk-i-  The  first  advantage  is  that  use  is  made 
of  prior  work,  saving  computation.  Second,  and  more  important,  is  that  Xk  involves  more  and 
more  data,  hence  more  and  more  storage  as  k  increases.  The  storage  of  Lk  requires  only  N2  real 
numbers. 

Since  the  interference  statistics  are  expected  to  change  with  time,  we  would  not  want  to  use 
the  scheme  described  above  without  modification.  But  interference  statistics  change  slowly  in 
comparison  with  the  duration  of  a  block  (or  else  we  would  not  be  able  to  hold  the  weight  vector 
constant  over  that  duration);  therefore  we  would  like  to  make  some  use  of  previous  blocks.  An 
often  used  compromise  is  to  weight  the  data  by  an  exponentially  decaying  window.  An  easy  way  to 
accomplish  this,  using  a  “forgetting  factor”  of  a,  is  to  use  aL*_i  in  place  of  Lk- 1  in  the  previous 
algorithm  for  updating  Lk . 

Since  L  can  be  updated  recursively,  we  can  now  update  L  for  each  new  vector  sample  as  it 
becomes  available  3  and  preserve  a  snapshot  of  L  =  Lk  after  the  update  that  represents  the  end 
of  the  kth  block.  This  makes  the  update  rule  very  simple.  Let  A  be  the  new  vector  sample  and 
let  L0id  be  the  Cholesky  matrix  for  all  previous  data.  Let  a  be  the  forgetting  factor.4  A  unitary 
matrix  Q  of  size  (N  +  1)  x  (N  - hi)  is  used  in  the  update: 

[aLold\X]Q^[L 

netulQ]  • 

Now  the  mechanization  of  the  unitary  transformation  can  be  discussed.  As  previously  ex¬ 
plained,  we  can  postmultiply  by  a  sequence  of  simple  unitary  transformations  to  accomplish  the 
same  effect  as  a  single  more  elaborate  transformation.  Since  the  purpose  is  to  introduce  zeroes  into 

3  In  an  earlier  study  by  Steinhardt  [5],  we  examined  the  approach  of  updating  L  for  many 
new  samples  at  once.  The  appropriate  unitary  transformation  for  multiple  sample  updates  is  a 
Householder  transformation.  Householder  transformations  require  less  arithmetic  than  the  ap¬ 
proach  taken  in  this  report,  but  they  have  not  yet  been  shown  to  have  a  simple  realization  in  highly 
parallel  systolic  hardware. 

4  a  is  now  chosen  based  on  updating  L  on  a  per-sample  basis  rather  than  on  a  per-block  basis. 
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the  tacked-on  column  X  of  the  matrix  [aL0/<f |X]>  we  will  use  unitary  transformations  which  each 
zero  out  one  additional  element. 

Two  types  of  unitary  transformations  are  used.  One  type  makes  the  leading  nonzero  element 
of  the  tacked-on  vector,  which  is  probably  a  complex  number,  become  a  real  number.  The  unitary 
matrix  Qe  which  accomplishes  this  is  of  the  form 


Qe  = 


‘  1 

0 

0  . 

.  0 

0 

0 

1 

0  . 

.  0 

0 

0 

0 

1  . 

.  0 

0 

0 

0 

0  . 

.  1 

0 

0 

0 

0  . 

.  0 

e* 

that  is  to  say,  Qe  is  an  (iV+ 1)  x  (N+ 1)  identity  matrix  with  the  lower  right  diagonal  element  replaced 
by  e*Q .  Postmultiplying  [aL0jd|X]  by  Q$  leaves  most  of  the  matrix  unchanged  but  multiplies  the 
last  column  (the  tacked-on  vector)  by  eJ0.  We  must  choose  0  such  that  multiplying  the  leading 
element  of  the  tacked-on  vector  by  e i6  makes  the  product  real. 

The  other  type  of  unitary  transformation  used  is  called  a  Givens  transformation  and  is  denoted 
by  Q#.  The  form  of  the  matrix  Q $  is  derived  from  an  (iV-f  1)  x  (iV  +  1)  identity  matrix  by  changing 
four  elements.  The  mth  and  (N  +  l)-st  diagonal  elements  become  cos<£.  The  element  in  row  (7V  + 1) 
and  column  m  becomes  sin0.  The  element  in  column  (JV  4-  1)  and  row  m  becomes  -  sin^>.  This  is 
illustrated  here  for  m  =  2: 


Q<t>  = 


'  1 

0 

0 

...  0 

0 

0 

cos  o 

0 

...  0 

—  sin  <j> 

0 

0 

1 

...  0 

0 

0 

0 

0 

...  1 

0 

0 

sin<£ 

0 

...  0 

COS0 

When  we  postmultiply  by  Q#  only  the  mth  column  and  the  tacked-on  ( N  +  l)-th  column  are 
affected.  Suppose  that  the  elements  in  the  same  row  in  those  two  columns  are  l  and  x.  After  the 
transformation 

l '  <=  l  cos  (j)  4-  x  sin  0,  and 

x'  <=  x  cos  (j)  —  l  sin  (j). 

If  l  and  x  are  real,  this  may  be  viewed  as  a  geometric  rotation  through  the  angle  (j)  in  the  plane 
defined  by  l  and  x  axes. 

By  an  alternation  of  postmultiplications  by  unitary  matrices  of  these  two  types,  the  entire 
tacked-on  column  can  be  zeroed  out.  This  progression  is  illustrated  in  Figure  2-2.  The  symbolism 
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in  this  figure  uses  ordinary  letters  to  indicate  complex  quantities  and  letters  with  an  overbar,  a,  to 
indicate  quantities  that  must  be  real.  In  addition,  if  a  column  in  a  matrix  on  the  left  is  in  boldface, 
that  column  is  affected  by  the  transformation.  Other  columns  are  not  affected.  An  example  with 
N  =  4  is  illustrated. 

When  the  mechanics  of  the  two  types  of  unitary  transformation  is  compared,  there  is  a  surpris¬ 
ing  commonality.  For  the  Qq  we  look  at  the  leading  nonzero  element  of  the  tacked-on  column  and 
determine  the  “angle”  9  that  will  make  this  leading  element  real.  Because  the  column  is  composed 
of  complex  numbers,  the  determination  of  9  depends  on  two  real  numbers  that  are,  respectively, 
the  real  part  and  the  imaginary  part  of  the  leading  element  x*: 


9  —  arctan 


-Im(xj) 

Re(xi) 


When  cos  6  and  sin  6  are  available,  the  transformation  Qq  is  applied  to  the  remainder  of  the  tacked- 
on  column,  which  we  may  consider  as  two  columns  of  real  numbers  that  are  its  real  and  imaginary 
parts.  In  terms  of  the  two  columns,  the  ith  Qq  transformation  takes  the  form: 


0  0 


0  0 
Re(xi)  Im(xi) 

Re(xi+ 1)  Im(xi+\) 


cos  9  —  sin  9 
sin  9  cos  9 


Re(xjv)  Jm(xjv) 


On  the  other  hand,  for  the  zth  Givens  transformation  Q^,  the  determination  of  the  “angle”  <j> 
is  made  on  the  basis  of  the  leading  elements  of  two  columns.  is  the  leading  element  of  column 
i  and  xt-  is  the  leading  element  of  the  tacked-on  column.  The  leader  of  the  tacked-on  column  is 
real  as  a  result  of  the  preceding  Qq  transformation.  Now  assume  that  is  real  -  therefore  by 
mathematical  induction  it  will  remain  real  because  all  the  operations  involved  in  updating  it  will 
produce  real  numbers.  <j>  is  determined  by 


<j>  ~  arctan 


which  is  exactly  analogous  to  the  determination  of  9  for  the  Qq  transformation.  Then,  when  cos  <p 
and  sin  <j>  are  available,  the  Givens  transformation  is  applied  to  the  only  two  columns  affected  by 
it.  The  Givens  transformation  may  be  written  in  terms  of  these  two  columns: 
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Figure  2-2.  Eight  successive  unitary  transformations  which  zero  out  the  four  element 
tacked-on  vector. 
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0 


0 


0  0 

X{ 

^i+l,i  ^t+1 

iNj  xn 


cos  (f>  —  sin  <\> 

sin  <£  cos  0 


This  shows  that  the  mechanization  of  the  Givens  transformation  is  almost  exactly  the  same  as 
the  mechanization  of  the  Qq  transformation.  There  is  one  difference.  After  dealing  with  the  leaders, 
the  remainder  of  the  Givens  transformation  is  applied  to  two  columns  of  complex  numbers,  whereas 
the  Q$  transformation  is  applied  to  two  columns  of  real  numbers.  However,  we  can  consider  the 
Givens  transformation  applying  separately  to  the  real  parts  of  its  two  affected  columns  and  to  the 
imaginary  parts  of  its  two  affected  columns.  In  these  terms  a  Givens  transformation  is  identical 
to  two  Qe  transformations.  The  unitary  transformations  needed  to  zero  out  a  tacked-on  column 
appended  to  a  (triangular)  Cholesky  matrix  are  all  identical  and  can  be  implemented  using  identical 
hardware.  In  the  subsequent  sections  a  design  for  such  hardware  in  detail  will  be  described. 

It  is  instructive  to  estimate  the  total  workload  involved  in  using  the  algorithm  we  have  described 
to  zero  out  a  single  64-element  vector.  Each  element  in  the  Cholesky  matrix  is  modified  using  three 
rotations.  To  carry  out  each  rotation  using  a  conventional  computer  architecture  would  appear 
to  require  two  load  instructions,  four  multiplications,  an  addition,  a  subtraction,  and  two  stores, 
about  ten  instructions  per  rotation.5  Therefore  3 (N  +  (N  —  1)  +  (N  —  2)  +  . . .  +  1)  =  %N(N  +  1)= 
6240  rotations,  equivalent  to  about  62,000  instructions,  are  involved  per  new  sample. 

Several  authors  have  proposed  parallel  computation  to  compute  the  update  of  a  Cholesky 
matrix.  The  best  known  architecture  is  a  systolic  array  of  computing  elements  arranged  in  a 
triangular  mesh  [4].  For  a  given  number  of  degrees  of  freedom,  jV,  such  a  triangular  array  uses 
N(N  —  l)/2  processors.  This  limits  practical  application  of  the  triangular  array  architecture  to 
very  small  N. 


5  The  steps  involved  to  rotate  the  leaders  are  different  and  more  involved,  but  we  are  counting 
them  as  if  they  were  the  same. 
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3.  CORDIC  REALIZATION  OF  GIVENS  TRANSFORMATION 


In  this  section,  an  algorithm  is  described  for  coordinate  rotation  that  is  well  suited  for  digital 
realization,  known  as  CORDIC.  The  acronym  stands  for  coordinate  rotation  digital  computation. 
The  basic  idea  was  first  published  in  the  1950s  by  Voider  [6]. 

Suppose  a  vector  points  from  the  origin  to  a  point  whose  coordinates  are  (x,t/),  and  that  this 
vector  is  to  be  rotated  to  a  new  point  with  coordinates  (x',yf),  such  that  the  angle  between  the 
new  and  old  points  is  £.  \Ve  have 

x '  =  (cos£)(x  —  y  tan£),and 

y'  =  (cos£)(y  +  xtan£)- 

This  involves  four  multiplications.  However,  there  are  many  special  angles  for  which  some  of  the 
multiplications  would  simplify  to  shifts.  We  will  concentrate  on  the  special  angles  for  which 

tan^  =  ±2_l/. 

Then  the  multiplications  by  tan  £„  become  right-shifts  by  v  bit  positions.  For  fixed  i/,  the  two  special 
angles  are  of  the  same  magnitude  but  opposite  sign  and  therefore  they  have  the  same  cosine. 

The  first  step  of  the  CORDIC  algorithm  is  to  realize  any  arbitrary  angle  £  by  a  sequence  of 
rotations  either  forward  or  backward  by  v  =  0, 1, . . . ,  oo.  Let  pv  =  ±1  determine  whether  a 

particular  “minirotation”  is  forward  or  backward.  Thus, 

oo 

£  = 

i/=0 


and  the  rotation  of  (x,  y)  through  the  angle  £  is  accomplished  by  the  sequence  of  steps  in  Figure  3-1. 

The  key  step  in  development  of  the  CORDIC  approach  is  to  recognize  that  the  multiplications 
by  cos^j/  may  all  be  collected  together  into  a  single  constant 

oo 

K  =  JJ  cos 

i/=0 


that  is  independent  of  the  overall  angle  £  by  which  we  rotate.  Thus  Figure  3-1  can  be  revised  to 
the  form  of  Figure  3-2. 

The  CORDIC  algorithm  is  composed  of  “stages” — most  of  the  stages  perform  minirotations  and 
a  “last”  stage  performs  a  gain  correction.  For  any  angle,  all  the  minirotations  are  employed,  but 
each  is  employed  in  either  a  clockwise  or  counterclockwise  direction.  Although  an  infinite  number 
of  minirotations  are  called  for  in  an  exact  realization  of  the  CORDIC  algorithm,  it  is  practical  to 
use  a  finite  number  of  stages,  in  the  same  sense  that  it  is  practical  to  make  a  binary  adder  or 
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x-2  °p0y  \ 
y  +  2~°p0x  J 
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<=  COS  f  i 


x-2  lpiy  \ 
y  +  2~lpxx  J 
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<=  COS  f  2 


x-2  2p2y  \ 

y  +  2~2p2x  J 


Figure  3-1.  Achieving  rotation  through  £  using  minirotations  through  the  special 
angles  ±£„. 


Figure  3-2.  Rotation  through  f  using  minirotations  through  the  special  angles 
using  one  correction  gain. 


16 


multiplier  with  a  finite  number  of  stages —  the  higher  numbered  stages  contribute  very  little  to  the 
accuracy  of  the  angle  specification.  The  correction  stage  is  a  multiplication  by  the  fixed  quantity 
K,  not  a  general- purpose  multiplier.  The  exact  value  of  K  depends  on  how  many  CORDIC  stages  are 
used.  As  Table  3-1  shows,  K(umax)  converges  quite  quickly.  (In  the  MUSE  application,  the  fixed 
multiplication  by  K<  designed  into  each  CORDIC  circuit,  is  combined  with  another  multiplication 
called  the  forgetting  factor.) 


TABLE  3-1. 

Correction  Constant  Required  when  the  Last  Stage  Rotates  by  arctan2”‘/mai 


Vmax 

K 

Vmax 

K 

0 

0.707106781 

8 

0.607254479 

l 

0.632455532 

9 

0.607253321 

2 

0.613571991 

10 

0.607253031 

3 

0.608833912 

11 

0.607252959 

4 

0.607648256 

12 

0.607252941 

5 

0.607351770 

13 

0.607252936 

6 

0.607277644 

14 

0.607252935 

7 

0.607259112 

15 

0.607252935 

A  CORDIC  method  therefore  achieves  rotation  without  using  any  of  the  trigonometric  functions 
and  without  explicit  multiplications.  If  the  angle  by  which  we  wish  to  rotate  the  pair  (x,  y)  is  known 
in  advance,  we  can  determine  the  set  of  controls 

(Pt  >  t  =  0?  •  •  •  >  Vmax  ) 

that  are  each  represented  by  a  single  bit.  This  may  be  thought  of  as  representing  the  angle  &  using 
digits  pi  in  an  unconventional  number  system,  a  number  system  different  from  binary,  decimal,  or 
any  of  the  radix  systems  that  are  commonly  used.  Each  pi  can  be  stored  in  a  flip-flop  in  the  stage 
whose  direction  of  rotation  it  is  to  control. 

However,  in  the  Givens  transformation  application  the  angle  of  rotation  is  not  known  in  ad¬ 
vance.  A  pair  (x,  y)  is  given  and  we  must  rotate  that  pair  through  the  angle  such  that  the  resulting 
rotated  pair  takes  the  form  (x',0).  This  operation  is  called  vectoring.  Then  we  must  rotate  some 
number  of  other  pairs  through  the  same  angle.  There  is  no  need  to  know  what  the  angle  is,  so  long 
as  we  are  able  to  rotate  by  that  angle.  Therefore  what  we  really  need  for  vectoring  is  an  algorithm 
to  determine  the  CORDIC  controls  pi.  A  major  advantage  of  the  CORDIC  algorithm  is  that  the  same 
circuit  that  is  used  for  rotating  may  be  used  for  vectoring. 

Consider  just  the  first  CORDIC  stage,  for  which  the  special  angle  is  either  45°  or  —45°.  Since 
the  purpose  is  to  rotate  the  input  (x,  y)  toward  the  x  axis,  if  y  is  above  the  axis  then  we  should 
rotate  “down,”  and  if  y  is  below  the  axis  then  we  should  rotate  “up.”  Therefore 
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Po  =  sgn(x)sgn(y). 


Once  po  is  determined  we  compute  the  effect  of  the  first  stage  on  x  and  y  and  pass  (x,  y)  as  modified 
to  the  second  stage.  Here  again  the  rule  is  to  rotate  down  if  y  is  above  the  axis  and  up  if  y  is  below 
the  axis. 


pi  =  sgn(x)sgn(y). 

In  this  way,  the  determination  of  the  controls  is  quite  simple. 

Pi  =  sgn(x)sgn{y),  i  =  0, . . . ,  vmax. 

These  controls,  once  determined,  are  saved  in  the  flip-flops  of  the  specialized  stages  and  are  used  to 
control  those  stages  for  the  succeeding  (x,  y)  pairs  that  are  to  be  rotated  through  the  same  angle. 

In  Figure  3-3  the  concept  of  a  CORDIC  circuit  made  up  of  independent  stages  is  shown.1  The 
inputs  (x,  y)  are  modified  by  minirotations  as  they  proceed  from  stage  to  stage.  This  circuit  has 
the  virtue  of  natural  pipelining.  If  registers  are  placed  between  the  stages  where  there  are  dashed 
lines  in  Figure  3-3,  then  a  new  rotation  problem,  involving  a  new  pair  ( x,y ),  can  be  started  by 


the  circuit  as  soon  as  the  preceding  pair  has  been  latched  at  the  output  of  the  first  stage;  then  yet 
another  rotation  may  be  started  when  the  first  two  pairs  have  been  latched  at  the  output  of  the' 
second  and  first  stages,  respectively,  etc.  Note  that  rotation  may  follow  vectoring  in  this  pipelined 

1  In  Figure  3-3  the  controls  pi  are  shown  as  if  they  depended  only  on  the  sgn(y).  The  simplified 
figure  is  meant  only  to  illustrate  the  pipelining  aspect  of  CORDIC  circuits.  In  fact  the  controls 
depend  on  sgn(y)sgn(x). 
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fashion.  Only  an  addition  and  a  subtraction  need  to  be  performed  in  each  stage.  These  operations 
can  be  carried  out  very  rapidly  in  digital  logic.  Therefore  the  rate  at  which  a  CORDIC  circuit  can 
begin  new  rotation  problems  is  quite  high.  By  contrast  the  time  required  to  complete  any  given 
rotation  is  proportional  to  the  number  of  CORDIC  stages  provided. 
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4.  CORDICS  IN  SYSTOLIC  ARRAY  FOR  UPDATING  A  CHOLESKY 

MATRIX 


The  pipelined  CORDIC  circuit  is  ideally  used  as  one  of  a  large  number  of  computing  elements 
operating  in  parallel  in  a  systolic  array.  The  term  “systolic”  originated  in  human  biology  to  describe 
the  system  of  circulation  of  the  blood  through  blood  vessels  and  various  organs.  In  a  systolic 
computing  system  the  data  is  “pumped”  from  computing  organ  to  computing  organ  and  processed 
as  it  moves.  Ideally  in  a  systolic  system,  data  should  be  processed  as  soon  as  it  becomes  available, 
so  it  must  be  available  at  the  right  place  at  the  right  time.  This  means  that  the  algorithmic  and 
the  architectural  considerations  are  very  closely  intertwined. 


4.1  Supercell 

In  this  section  we  will  show  how  a  large  number  of  the  pipelined  CORDIC  circuits  can  be  used 
in  a  parallel  systolic  pipelined  system  to  update  a  Cholesky  matrix.  One  CORDIC  circuit  will  be 
assigned  to  perform  each  of  the  N  Q$  transformations,  and  two  CORDIC  circuits  will  be  assigned  to 
perform  each  of  the  N  Q<p  transformations.  Therefore  three  CORDIC  circuits  are  needed  to  perform 
all  the  operations  associated  with  updating  one  of  the  N  columns  of  the  Cholesky  matrix  and  3 N 
CORDIC  circuits  are  needed  in  all. 

A  common  data  format  will  be  used  throughout.  Data  that  are  vectors  of  complex  numbers 
are  presented  to  circuitry  sequentially,  and  the  time  during  which  one  complex  word  enters  a 
CORDIC  circuit  will  be  called  a  microcycle.  Every  vector  that  enters  a  circuit  will  have  a  leader,  its 
first  element,  and  some  number  of  followers,  all  the  other  elements.  The  CORDIC  circuits  perform 
vectoring  on  the  leader  and  rotation  through  the  same  angle  for  the  followers.  A  vector  composed 
of  K  elements,  for  example,  will  flow  into  a  CORDIC  circuit  during  K  consecutive  microcycles.  The 
elements  making  up  such  a  vector  will  flow  out  of  the  CORDIC  circuit  at  the  same  rate  that  they 
entered,  one  element  per  microcycle.  Although  the  elements  are  modified  by  passage  through  the 
CORDIC  they  retain  their  identity  and  their  order,  and  the  leader  on  input  becomes  the  leader  on 
output.  The  CORDIC  circuit  has  a  delay  due  to  pipelining,  which  is  called  its  latency. 

The  (?0-type  transformations  are  simply  the  passage  of  a  vector  through  a  CORDIC  circuit. 
The  action  of  the  CORDIC  circuit  on  the  leader  is  a  phase  change  such  that  the  leader  becomes  real. 
The  action  of  the  CORDIC  circuit  on  the  followers  is  to  change  their  phases  by  the  same  amount. 
A  CORDIC  circuit  used  in  this  manner  will  be  called  a  0-CORDIC. 

While  a  Qq  transformation  is  a  phase  change  applied  to  a  single  complex  vector,  a  Q $  trans¬ 
formation  is  a  rotation  of  pairs  of  complex  numbers  through  a  real  angle.  To  accomplish  a  Q $ 
transformation  we  use  two  CORDIC  circuits,  which  are  called  the  master  0- CORDIC  and  the  slave 
0- CORDIC,  respectively.  The  master  0-CORDIC  deals  with  real  parts  of  complex  data,  while  the 
slave  0-CORDIC  deals  with  corresponding  imaginary  parts.  Data  involved  in  a  Q $  transformation 
are  presented  to  the  two  CORDlCs  as  “vectors”  of  “complex”  elements  so  that  the  rotation  may 
be  accomplished  is  if  it  were  a  phase  change;  both  CORDlCs  rotate  using  the  same  angle.  These 
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vectors  are  presented  sequentially,  one  element  per  microcycle,  just  as  for  a  0-CORDl‘c.  Refer  to 
Figure  4-1.  The  input  vectors  for  0-CORDiCs  have  complex  elements  made  up  partly  from  the 
output  of  a  0-CORDIC  and  partly  from  one  column  of  the  Cholesky  matrix.  The  master  </>-CORDlC 
gets  the  real  part  of  its  input  vector  from  the  real  part  of  the  Cholesky  matrix  column  and  gets 
its  imaginary  part  from  the  real  part  of  the  0-CORDIC  output.  The  slave  <jh CORDIC  gets  the  real 
part  of  most  of  its  input  vector  from  the  imaginary  part  of  the  Cholesky  matrix  column  and  gets 
the  corresponding  imaginary  parts  from  the  imaginary  part  of  the  0-CORDIC  output.  However,  the 
leader  of  the  vector  that  is  input  to  a  slave  <jf>-CORDlC  is  a  special  case,  which  is  discussed  below. 
The  outputs  of  the  two  types  of  0-CORDlCs  are  then  reassembled.  The  real  parts  of  the  <^CORDlC 
outputs  become  the  real  and  imaginary  parts,  respectively,  of  the  updated  Cholesky  matrix  column. 
The  imaginary  parts  of  the  </>-CORDlC  outputs  become  the  real  and  imaginary  parts,  respectively, 
of  the  updated  tacked-on  vector.  A  more  complete  explanation  follows. 

Note  that  CORDIC  circuits  are  used  for  two  quite  different  meanings  of  “rotation,”  one  being 
the  phase  modification  of  a  complex  number  and  the  other  being  the  coordinate  rotation  of  a  two- 
element  vector.  In  the  latter  case,  although  the  two  coordinates  may  be  complex  numbers,  the 
angle  through  which  they  are  rotated  is  real,  and  the  same  real  rotation  angle  is  used  on  the  real 
and  imaginary  components,  respectively. 

First,  consider  the  zth  transformation  Q$.  Its  input  is  the  tacked-on  vector  with  components 
[xj,  X{+ 1, . . . ,  xjv],  all  complex  numbers.  These  are  fed  into  a  0-CORDIC  one  pair  at  a  time,  beginning 
with  the  pair  (i?e(xj),  Im(xi)).  The  first  pair  is  marked  as  a  “leader”  with  a  special  indicator  that 
causes  the  CORDIC  circuit  to  determine  and  set  the  pi  controls  as  it  passes  through  the  successive 
rotation  stages  (vectoring).  Thus  the  modified  Xi  that  emerges  from  the  0-CORDIC  is  rotated  in 
phase  so  that  it  is  real — its  imaginary  component  is  zeroed.  The  “followers”  X{+ i,...,xjv  are 
rotated  by  the  same  phase  as  the  leader,  but  they  generally  emerge  from  the  0-CORDIC  as  complex 
numbers. 

It  is  helpful  to  think  of  the  tacked-on  vector  as  being  pumped  through  the  0-CORDIC  one 
element  at  a  time,  and  to  think  of  the  CORDIC  as  a  pipelike  organ.  As  many  elements  as  the 
CORDIC  has  stages  will  be  inside  the  pipe  at  any  time. 

In  the  ith  transformation  Q#  the  data  being  pumped  through  the  transformation  are  two 
columns,  [Z^,  Zi+i,i, . . . ,  Z^i]  and  [x*,  x*+i, . . . ,  xjv]-  These  are  composed  of  complex  numbers  except 
that  the  pair  (lij,Xi)  is  real.  Things  are  arranged  so  that  Z*.,*  and  x are  available  simultaneously. 
One  CORDIC,  the  master  (^-CORDIC,1  gets  its  two  inputs  from  the  real  parts  of  the  two  input  streams. 
The  slave  (fh CORDIC  takes  its  two  inputs  (with  one  exception)  from  the  imaginary  parts  of  the  two 
input  streams.  The  pair  (Z^x*)  is  marked  as  a  leader  so  that  as  it  passes  through  the  master 
(jy CORDIC  it  determines  the  angle  controls  pi  for  that  CORDIC,  which  are  to  be  used  to  rotate  all 
the  follower  pairs  (Re(lk}i),  Re(xk))  which  pass  through  the  master  </>-CORDlC.  The  same  angle 
controls  are  to  be  used  in  the  slave  <£-CORDlC  to  rotate  the  pairs  (Im(lk%i),Im(xk))-  The  slave 
^-CORDIC  is  not  needed  for  (/ra(Zj}j),  Im(xi))  because  this  pair  is  (0,0).  Therefore  it  is  “tricked” 
into  setting  its  angle  controls  identical  to  those  of  the  master  0-CORDIC  by  giving  it  the  leader 

1  The  terms  master  and  slave  are  left  over  from  an  earlier  usage  and  are  not  meaningful. 
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pair  (Re(liti),  Re(xi))  in  place  of  (Im(liti),  7m(x,)).  The  outputs  of  the  (^-CORDlCs  are  an  updated 
column  of  the  Cholesky  matrix  and  an  updated  tacked-on  vector  with  its  ith  element  zeroed. 

In  this  way,  the  use  of  three  CORDIC  circuits  accomplishes  the  pair  of  transformations  needed 
to  zero  out  one  element  in  the  tacked-on  vector.  This  three-CORDIC  configuration  is  called  a 
supercell.  In  Figure  4-1  we  illustrate  how  the  0-CORDIC  and  the  master  and  slave  0-CORDlCs  are 
interconnected  and  how  the  pipelining  of  each  individual  CORDIC  circuit  extends  naturally  to  the 
pipelining  of  the  supercell. 

144246-3 

INPUT  TACKED-ON 


OUTPUT  • 
TACKED-ON  COLUMN 


Figure  4-1.  A  CORDIC  supercell. 

Figure  4-1  shows  how  the  column  [lij,  U+i,i, . . . ,  l^.i]  can  be  stored  within  the  supercell  in  local 
memory.  On  the  other  hand  [xj+i, . . . ,  xyvj  must  be  passed  to  another  supercell  where  it  is  needed 
by  the  (i  +  l)st  Qe  transformation  and  subsequent  Q$  transformation.  Since  the  leader  of  the 
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tacked-on  vector  is  zeroed  out  by  passing  through  a  supercell,  the  element  immediately  following 
it  in  the  sequence  becomes  the  new  leader  element. 


4.2  Simple  Systolic  Array 

Now  that  there  is  a  configuration  for  a  supercell  that  can  update  one  column  of  a  Cholesky 
matrix,  a  complete  systolic  array  using  N  such  supercells  can  be  configured.  Figure  4-2  shows  such 
a  systolic  array,  except  that  the  local  memory  in  each  supercell,  which  holds  the  column  that  the 
supercell  must  update,  is  shown  separately.  In  this  array  the  data  passed  from  supercell  to  supercell 
are  components  of  the  tacked-on  vector,  one  complex  word  at  a  time.  Thus,  at  the  input  to  the 
first  supercell  the  data  seen,  in  order  of  arrival,  are: 

.  .  .  |,  |xi,X2,X3,  .  .  .  ,Xat-i,Xat|,  |xi,X2,  .  .  . 

Each  block  consists  of  one  vector  sample  of  observed  interference  X_.  Note  that  as  soon  as  one 
block's  vector  has  been  pumped  into  the  first  supercell ,  it  is  ready  to  receive  another  block's  vector . 
The  time  interval  between  the  first  elements  of  two  successive  vectors  will  be  called  a  macrocycle. 

The  output  of  the  first  supercell  is  a  modified  tacked-on  vector.  Its  other  product,  the  updated 
first  column  of  the  Cholesky  matrix,  is  retained  within  the  supercell,  to  be  used  as  original  data  for 
the  next  update.  The  modified  tacked-on  vector  has  one  fewer  component  than  the  original.  Thus 
the  output  of  the  first  supercell  (input  to  the  second  supercell)  takes  the  form 

.  .  .  j ,  |6,  X2 ,  X3 , . . .  ,  Xyy _  1 ,  XN  j ,  1 6,  X2 ,  •  •  • 

where  the  b’s  are  blanks,  e.g.,  time  intervals  when  no  data  is  present.  The  second  supercell’s  output 
(the  third  supercell’s  input)  has  two  blanks  per  block: 

. . .  | ,  (6, 6,  X3, . . . ,  x^v_i ,  xjv|,  | h,b,  X3, . . . 

Strictly  speaking,  the  blanks  represent  zero-valued  samples.  But  these  are  not  values  that  just 
happen  to  be  zero —  they  are  inevitably  zero.  Thus  they  need  not  be  explicitly  sent  from  supercell  to 
supercell,  nor  involved  in  computation.  Therefore  the  first  supercell  is  busy  all  the  time,  the  second 
supercell  is  idle  for  one  microcycle  out  of  every  A,  the  third  supercell  is  idle  for  two  microcycles 
out  of  every  N,  and  the  last  supercell  is  idle  for  all  but  one  microcycle  out  of  every  N. 

4.3  Latency 

Since  a  CORDIC  circuit  is  pipelined,  there  will  be  a  circuit  latency.  This  means  that  there  is 
a  delay  between  presenting  an  element  at  the  input  and  getting  the  corresponding  result  at  the 
output.  Latency  is  essentially  the  amount  of  time  the  data  stays  in  the  “pipe.”  A  supercell  has  the 
latency  of  two  CORDIC  circuits.  We  will  use  r  to  represent  the  latency  of  a  supercell  in  microcycles. 


24 


144246-4 


INPUT 
N  WORDS 
PER  SAMPLE 


SUPER 

CELL 


MEMORY  FOR  COLUMN  N 


Figure  4-2.  Systolic  array  to  update  a  Cholesky  matrix. 


Thus  as  a  simple  example,  if  r  =  3  and  N  =  5  then  there  would  be  the  following  time  relationship 
among  the  five  (numbered  0  to  4)  supercells: 


-  0 
0-  1 
1  -  2 
2-3 


X\  3^2  ^3  X4  X5 

b  X2  £3  £4  £5 

b  b  £3  £4  £5 

b  b  b  X4  X5 


3-4: 


b  b  b  b  X5 


This  is  the  timing  relationship  for  successive  modifications  of  one  tacked-on  vector  as  it  is 
passed  from  supercell  to  supercell.  However,  because  each  supercell  begins  to  work  on  another 
update  as  soon  as  it  completes  the  previous  update,  a  much  busier  timing  relationship  would 
actually  be  observed: 
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This  illustrates  one  way  that  the  systolic  array  is  inefficient.  The  higher  numbered  supercells 
have  many  microcycles  in  which  they  are  not  doing  anything.  In  addition,  the  local  memory  used  is 
different  in  each  0-CORDIC  since  the  columns  of  L  are  of  different  length.  The  cells  must  either  be 
provided  with  different  amounts  of  memory,  or,  if  they  are  all  identical,  some  of  them  must  waste 
much  of  their  provided  memory. 

Despite  the  inefficiency,  we  should  note  the  excellent  characteristics  of  this  architecture  for 
updating  a  Cholesky  matrix: 

•  The  CORDIC  circuits  are  highly  specialized  for  rotation,  hence  they  are  efficient. 

•  The  CORDIC  circuits  are  naturally  pipelined,  hence  they  may  be  clocked  at  high 
speed. 

•  Data  is  passed  only  locally  from  processor  to  nearby  processor,  hence  there  is 
no  need  for  a  global  bus. 

•  All  CORDIC  circuits  are  operated  in  synchrony,  using  the  same  clocks. 

•  Data  arrives  where  it  is  needed  just  exactly  when  it  is  needed,  hence  there  is 
no  need  for  buffer  storage. 


4.4  Forgetting  Factor 

In  practice,  we  want  to  update  the  Cholesky  matrix  as  if  the  data  were  weighted  by  an  expo¬ 
nential  time  window.  This  can  be  accomplished  if  each  <f>-C ORDIC  attenuates  its  L  output  slightly 
by  a  forgetting  factor  a  before  reusing  it  in  the  next  update.  Since  the  CORDIC  must  have  a  cor¬ 
rection  gain  K  as  well,  we  propose  to  realize  a  single  correction  gain  of  Ka.  This  choice  has  some 
consequences.  Because  there  is  some  flexibility  in  the  choice  of  a,  we  can  choose  it  so  that  we 
can  multiply  by  Ka  with  minimum  hardware.  We  were  able  to  use  only  three  adder  stages  in  the 
eventual  choice,  a  =  0.99867003.  Each  tacked-on  vector  from  n  updates  earlier  affects  the  Cholesky 
matrix  as  if  it  had  been  multiplied  by  an.  Specifically,  a320  is  approximately  0.6532. 

We  would  prefer  that  the  corrected  gain  of  0-CORDlCs  be  unity,  and  that  the  corrected  gain  for 
(fhC ORDlCs  also  be  unity  when  correcting  the  gain  for  the  tacked-on  vector.  However,  it  simplifies 
the  hardware  design  if  we  use  only  one  correction  gain.  The  attenuation  of  the  tacked-on  vector 
due  to  the  forgetting  factor,  as  the  tacked-on  vector  is  passed  through  128  CORDlCs  in  series,  is 
slight  enough  that  it  may  be  shown  to  make  no  difference  in  nulling  performance. 
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5.  AN  EFFICIENT  SYSTOLIC  ARRAY  WITH  LATENCY  CONTROLLED 

INTERLEAVING 


In  the  previous  section  a  systolic  array  for  updating  a  Cholesky  matrix  was  introduced.  In  this 
section  a  major  improvement  that  increases  the  efficiency  of  the  systolic  array  is  explained. 

To  begin,  each  supercell  was  assigned  the  responsibility  for  updating  not  one  but  two  columns 
of  the  Cholesky  matrix.  The  lengths  of  the  two  assigned  columns  will  always  add  up  to  N  +  l.1 
This  is  done  as  in  Table  5-1. 


TABLE  5-1. 


Pairing  of  Columns  in  the  Same  Supercell 


Supercell 

Column  No. 

Length 

Column  No. 

Length 

0 

1 

N 

N 

1 

1 

2 

N-l 

N-l 

2 

2 

3 

N-2 

N-2 

3 

k 

k+1 

N-k 

N-k 

k+1 

N/2-1 

N/2 

N/2+1 

N/2+1 

N/2 

A  modification  of  the  simple  systolic  array  to  use  the  pairing  leads  us  to  Figure  5-1.  By  this 
pairing  we  obviously  equalize  the  workload  of  the  supercells,  since  each  supercell  is  required  to 
deal  with  two  leaders  and  N  —  1  followers  for  every  new  vector  sample  of  interference  that  comes 
along.  Only  half  as  many  supercells  are  needed  as  before.  Furthermore,  the  memory,  which  is  used 
in  0-CORDICs  to  store  the  columns  of  the  Cholesky  matrix  to  be  updated  in  that  supercell,  is  the 
same  size  in  all  the  supercells,  N  +  1  words  per  <£-CORDlC. 

This  modification  allows  us  to  retain  all  the  advantages  of  a  systolic  array.  There  are  two 
aspects  of  this  statement  that  need  clarification.  First,  because  the  tacked-on  vector  in  the  originally 
proposed  systolic  array  is  passed  from  supercell  k— 1  to  supercell  fc,  and  since  in  the  modified  systolic 
array  supercell  k  has  the  responsibilities  of  supercell  N  —  1  —  A:,  each  supercell  must  pass  data  both 
forward  and  backward.  For  example,  supercell  2  must  pass  its  output  to  supercell  3  but,  later, 
when  it  is  serving  the  role  of  supercell  N  —  3  it  must  pass  its  output  to  supercell  N  —  2  which  is 
really  supercell  1.  Yet  these  communications  are  still  only  local.  There  is  still  no  need  for  a  global 
bus. 

1  We  are  assuming  that  N  is  even.  A  similar  line  of  development,  not  taken  in  this  report,  can 
be  worked  out  if  N  is  odd.  We  would  use  the  first  supercell  for  the  longest  column  and  pair  the 
remaining  columns  so  that  their  lengths  added  up  to  N. 
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Figure  5-1.  Global  architecture  of  folded  systolic  array. 
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The  second  consideration  is  more  difficult.  In  a  systolic  array  we  would  like  the  data  to  get  to 
where  it  is  needed  exactly  when  it  is  needed.  It  would  not  he  suitable  if  supercell  k  were  passed 
tacked-on  vectors  from  supercell  k  +  1  and  also  from  supercell  k  —  1  during  the  same  microcycles. 
To  explore  this  issue  in  more  detail  we  need  to  take  account  of  the  supercell  latency  r  as  well  as 
the  workload  of  each  supercell.  Further,  we  must  consider  more  than  one  update  cycle  at  a  time. 
For  example,  supercell  0  would  devote  TV  microcycles  to  the  update  of  column  1  for  the  current 
tacked-on  vector,  then  devote  one  microcycle  to  the  update  of  column  TV  for  a  much  older  input.  In 
supercell  0,  suppose  that  the  current  tacked-on  vector  is  input  during  microcycles  0, 1, . . . ,  iV  —  1. 
When  supercell  0  serves  as  supercell  -AT  —  1,  it  works  on  a  tacked-on  vector  of  length  1.  Therefore 
we  must  leave  at  least  one  blank  microcycle  per  cycle.  It  would  be  best  if  the  one-element  vector 
arrived  exactly  as  the  supercell  becomes  available  during  microcycle  TV.  Then  the  whole  timing 
pattern  can  repeat  periodically  beginning  in  microcycle  TV  +  1  with  a  new  TV-element  tacked-on 
vector  during  microcycles  TV  +  1,  TV  +  2, . . . ,  2TV  followed  by  a  blank  in  microcycle  2 TV  -f-  1  during 
which,  as  supercell  TV  —  1,  a  new  one-element  tacked-on  vector  is  input. 

Letting  the  latency  be  an  unknown  quantity  r,  we  can  deduce  when  inputs  arrive  at  supercell  1, 
then  at  supercell  2,  etc.,  until  a  pattern  is  obvious  (see  Table  5-2).  The  pattern  allows  us  to  write 
the  list  for  any  generic  supercell,  Pj.,  and  for  the  generic  supercell  P/v-i-fc  that  is  its  alter  ego.  The 
microcycles  when  inputs  are  presented  are  congruent  mod  TV  +  1  to  the  lists  in  Table  5-2. 


TABLE  5-2. 


Microcycles  when  Each  Supercell  Input  Is  Busy 


Processor 

When  Inputs  Are  Presented 

Inputs/cycle 

Po 

0 . N-l 

N 

Pi 

r+1 . r-f  N-l 

N-l 

P2 

2(t+1) . 2r+N-l 

N-2 

Pk 

k(r+l) . kr+N-1 

N-k 

PN-l-k 

(N-l-k)(r+l),. . .  f(N-l-k)r+N-l 

k+1 

The  key  point  is  that  r  must  be  chosen  so  that  all  these  inputs  are  presented  to  supercells 
during  all  these  microcycles  without  any  data  collisions.  It  is  completely  adequate  to  derive  a 
condition  on  r  from  the  generic  supercell  that  serves  as  Pk  and  The  only  legitimate  time 

for  supercell  Pk  to  change  its  role  to  supercell  P/v-i-fc  is  just  after  microcycle  kr  +  TV  —  1,  when 
the  latter  must  accept  the  arrival  of  data  from  the  other  direction,  beginning  in  a  microcycle  whose 
number  is  congruent  to  (TV  —  1  —  k)(r  +  1)  mod  TV  +  1. 

{kr  +  TV  —  1)  4- 1  =  (TV  —  1  —  k)(r  +  1)  mod  (TV  +  1) 
kr  —  1  =  (—2  —  k)(r  +  1)  mod  (TV  +  1) 
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k(2r  +  1)  =  — (2r  4-  1)  mod  (N  +  1) 


( k  4-  l)(2r  +  1)  =  0  mod  ( N  4-  1) 

At  this  point  if  the  same  r  is  to  work  for  all  k  then  we  must  have 
2r  +  1  =  0  mod  (N  4-  1). 

The  above  equation  is  a  congruence.  We  can  write  it  as  an  equality  with  an  unknown  multiplier 
m 


2r  4-  1  =  m(N  4-  1). 


The  smallest  allowable2  latency  would  be  when  m  =  1,  namely 


r  = 


N 
2  * 


The  idea  of  choosing  a  latency  that  allows  all  the  various  partially  zeroed  tacked-on  vectors 
of  different  lengths  to  interleave  with  one  another  without  collisions  is  called  latency-controlled 
interleaving.  It  allows  us  to  “fold”  the  linear  systolic  array,  which  was  only  about  50  percent 
efficient,  into  a  half-size  linear  array  that  is  100  percent  efficient,  without  paying  a  penalty  in 
complication  of  control. 

Earlier,  when  the  inputs  to  <£-CORDICs  was  discussed,  we  arranged  that  which  comes  from 
a  local  memory,  and  which  comes  from  a  0-CORDIC,  should  arrive  at  the  <?!>-CORDlC  input  during 
the  same  microcycle.  This  is  also  accomplished  by  latency-controlled  interleaving,  but  of  a  much 
simpler  sort.  The  delay  from  one  update  cycle  to  the  next  is  TV  4-  1  microcycles.  Therefore,  as  lkti 
is  updated,  the  updated  /*,i  must  be  delayed,  by  a  combination  of  the  latency  of  the  0-CORDIC  and 
some  extra  delay,  so  that  it  appears  at  the  input  again  N  +  1  microcycles  later.  The  latency  of  a 
0-CORDIC  is  approximately  The  extra  delay  needed  can  be  provided  by  a  small  memory,  which 
has  about  N  4-  1  —  \  words. 

In  Figures  5-1  and  5-2  the  folded  systolic  array  is  shown.  Figure  5-1  is  the  global  architecture 
and  Figure  5-2  shows  the  communication  for  a  typical  supercell  and  its  neighbors.  The  necessary 
control  of  this  systolic  array  has  two  aspects.  First,  the  data  passed  from  CORDIC  cell  to  CORDIC 
cell  must  be  marked  with  an  indicator  of  whether  it  is  a  leader  or  a  follower.  Second,  since  a 
0-CORDIC  cell  can  take  its  input  from  two  other  supercells,  forward  or  backward  in  the  chain,  some 
control  must  tell  it  which  to  choose. 

A  simple  rule  accomplishes  both  controls.  At  the  input  to  the  first  supercell  (supercell  0),  an 
input  bit  called  direction  is  provided.  When  direction  =  true  the  leading  supercell  takes  its  input 
from  the  given  JV-element  tacked-on  vector.  When  direction  changes  to  false  the  leading  supercell 

2  Because  N  is  even,  a  parity  argument  will  show  that  m  must  also  be  odd. 
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Figure  5-2.  Local  communication  within  folded  systolic  array. 

takes  its  next  input  from  supercell  1  (acting  as  supercell  N  —  2).  The  change  in  direction  from  false 
to  true  or  from  true  to  false  serves  as  a  mark  of  an  element  that  is  a  leader,  causing  the  CORDIC  to 
set  its  controls  as  that  element  passes  through  the  stages. 

The  direction  bit  is  passed  systolicly  along  with  the  tacked-on  vector  through  both  types  of 
CORDIC  cells  and  emerges  from  the  <£-CORDlCs  to  be  used  in  the  next  supercell.  However,  the 
leading  edge  must  be  delayed  by  one  additional  microcycle  relative  to  the  tacked-on  vector  so  that 
the  change  in  direction  is  attached  to  the  first  element  not  zeroed  out.  The  direction  bit  travels  only 
forward  along  the  chain.  Direction  bits  output  during  processing  of  data  being  passed  backward 
would  be  redundant  and  are  ignored.  In  this  way,  the  control  of  the  folded  array  is  actually  no  more 
complex  than  the  control  of  the  original  systolic  array,  and  each  supercell  generates  the  control  for 
the  supercells  that  follow  it  in  the  chain. 
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6.  SIMULATION  TO  SELECT  PARAMETERS 


We  want  MUSE  to  be  able  to  demonstrate  adaptive  nulling  with  S/N  improvement  of  at  least 
50  dB.  Now  there  are  many  factors  in  a  system  that  control  the  achievable  improvement  in  S/N 
due  to  adaptive  nulling.  These  include  the  number  and  distribution  of  jammers,  their  strength  in 
comparison  to  thermal  noise  in  the  system,  the  antenna  element  gain  patterns,  mismatch  in  the  front 
ends  and  in  all  other  circuitry  prior  to  digitization,  the  adequacy  of  the  statistical  representation 
of  interference  (principally,  how  many  samples  used  to  characterize  it)  and,  finally,  the  accuracy  of 
the  digital  computation.  MUSE’s  design  should  not  rely  on  special  knowledge  of  the  environment, 
the  antenna  elements,  or  front  ends.  The  point  of  view  is  taken  that  MUSE  should  be  capable  of 
adapting  to  artificial  random  inputs  generated  with  only  the  constraint  that  50  dB  of  nulling  is 
possible.  Even  here,  however,  there  is  a  caveat.  Since  it  is  relatively  easy  to  provide  artificial  data 
whose  dynamic  range  will  overwhelm  any  fixed  design,  we  must  require  that  the  artificial  data  sets 
include  some  bound  on  the  largest  jammer  to  thermal  noise  ratio,  and  indeed,  this  bound  should 
be  only  a  little  bit  more  than  the  desired  50  dB  S/N  improvement  sought. 

With  this  goal,  we  have  been  able  to  select  MUSE  parameters  by  extensive  numerical  simu¬ 
lation.  Each  hardware  configuration  under  consideration  was  simulated  exactly  at  the  “bit-level” 
with  the  chosen  number  of  stages,  word-length,  etc.  Simulated  data  sets  were  contrived  with  any 
selected  number  of  “jammers,”  such  that  exactly  50  dB  of  S/N  improvement  (or  any  other  desired 
amount)  was  theoretically  possible  for  that  data  set.  The  “condition  number”  of  the  contrived 
data  was  also  controlled  so  as  to  be  realistic.  Except  for  these  two  constraints,  the  data  sets  were 
random  and  several  data  sets  with  the  same  condition  number  and  theoretical  S/N  improvement 
potential  were  tried  for  each  experiment.  In  appendix  A,  the  mathematical  method  by  which  we 
can  produce  data  from  a  random  scenario  wTith  these  two  constraints  is  explained.  The  Cholesky 
matrix  determined  by  MUSE  for  each  data  set  was  used  to  determine  weights  and  the  weights  were 
used  to  determine  the  actual  S/N  improvement.  Therefore,  the  loss  in  performance  due  to  the  finite 
word-length,  finite  number  of  CORDIC  stages,  rounding  and  truncation,  etc.  was  determined. 

Based  on  these  simulations,  four  finite  word-length  parameters  could  be  determined:  the  word- 
length  of  internal  CORDIC  registers,  the  word-length  for  the  tacked-on  vector,  the  word-length  for 
the  Cholesky  matrix,  and  the  number  of  CORDIC  stages.  The  “forgetting  factor”  was  also  studied 
and  is  described  elsewhere. 

The  heart  of  the  simulation  is  a  subroutine  that  simulates  a  CORDIC  process  with  integer 
arithmetic.  Because  a  typical  CORDIC  stage  requires  scaling  data  by  2“l  and  adding  it  to  unsealed 
data,  the  word-length  required  for  a  precise  execution  of  the  equations  in  Figure  3-2  would  increase 
the  number  of  bits  to  the  right  of  the  binary  point  by  i  bits  in  the  zth  stage.  Thus  we  need  rounded 
arithmetic  in  almost  every  stage.  However,  such  rounding  introduces  computational  noise  in  every 
stage,  and  the  total  noise  could  be  excessive.  Guard  bits  introduced  at  the  low  end  of  each  word 
reduce  the  rounding  noise.  Experiments  showed  that  two  guard  bits  reduced  the  rounding  noise 
adequately. 
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It  was  also  necessary  to  consider  overflow.  Of  course,  overflow  in  a  CORDIC  might  be  inevitable. 
If  x  and  y  are  both  near  their  maximum  allowed  values  at  input  and  are  rotated  by  45°,  then  a 
component  of  the  result  can  exceed  the  maximum  by  a  factor  of  almost  \/2.  This  kind  of  overflow 
must  be  prevented  by  appropriate  scaling  of  the  inputs,  a  consideration  common  to  all  fixed-point 
signal  processing  projects.  However,  a  second  kind  of  overflow  is  also  a  concern — overflow  within 
an  interior  CORDIC  stage.  This  was  “prevented”  in  connection  with  realizing  the  combination  of  the 
CORDIC  gain  compensation  and  forgetting  factor.  The  required  fixed  multiplication  (1242/2048)  is 
achieved  by  using  three  shift-and-add  circuits,  e.g., 

1242  (27  +  23  +  21)(23  +  l) 

2048  “  211 

The  multiplication  is  split  between  a  factor  of  138/256  at  the  input  to  the  CORDIC  (the  guard  bits 
are  created  at  this  point)  and  9/8  at  the  output  of  the  CORDIC  (the  guard  bits  are  dropped  at  this 
point).  The  reason  for  splitting  the  factor  in  this  way  is  to  ensure  that  overflows  can  only  occur  in 
the  CORDIC  circuit’s  interior  in  cases  when  overflows  would  be  inevitable  due  to  simply  storing  the 
result  of  an  accurate  rotation  in  the  available  word-length. 

Because  the  simulation  of  the  CORDIC  takes  exact  account  of  the  finite  word-length  arithmetic 
process,  the  results  it  produces  are  bit-by-bit  identical  to  the  results  expected  from  actual  hardware. 

The  simulated  CORDIC  subroutine  is  called  with  each  pair  of  inputs,  and  with  a  flag  telling 
it  whether  this  is  a  leader  pair  for  which  the  controls  are  to  be  set  or  a  follower  pair  for  which 
previous  controls  are  to  be  used.  The  subroutine  is  called  three  times  within  another  subroutine 
that  simulates  a  supercell  for  a  pair  of  inputs.  An  L-Update  subroutine  calls  the  supercell  subroutine 
64  +  63  +  62  + . . .  +  2  +  1  =  2080  times  for  each  tacked-on  vector.  The  L-Update  subroutine  is  called 
once  for  each  vector  sample  in  the  data  set.  This  simulation  ignores  the  latency  and  other  timing 
issues,  since  they  have  no  bearing  on  the  relations  between  parameter  choices  and  performance. 

The  contrived  data  sets  comprised  exactly  N  =  64  vector  samples  for  which  the  ideal  Cholesky 
matrix  was  predetermined.  Ideal  weights  were  computed  using  floating  point  arithmetic  for  this 
ideal  Cholesky  matrix  to  confirm  that  the  predetermined  S/N  improvement  was  correct.  Then  these 
64-sample  data  sets  were  repeated  as  MUSE  inputs  several  times  in  succession,  and  snapshots  of 
L  were  taken  after  every  block  of  64  samples.  The  computed  Cholesky  matrices  were  then  used  to 
compute  MUSE  weights  using  accurate  arithmetic.  Any  loss  in  S/N  improvement  when  the  MUSE 
weights  are  used  in  place  of  the  ideal  weights  could  therefore  be  attributed  to  the  computational 
error  in  computing  L  with  finite  word-length  arithmetic. 

Another  method  of  computing  weights  will  be  described  in  section  8  in  which  the  linear  equa¬ 
tions  for  the  desired  weights  are  solved,  after  L  is  already  computed,  using  the  same  finite  word- 
length  CORDIC  arithmetic.  Therefore  nulling  based  on  weights  computed  that  way  should  generally 
show  an  additional  loss  in  S/N  improvement.  This  was  also  studied. 

In  Figure  6-1,  the  abscissa  gives  the  condition  number  of  a  possible  data  set  and  the  ordinate 
gives  the  S/N  improvement  possible  for  that  data  set.  The  solid  curve  gives  the  upper  bound  of 
the  set  of  possible  scenarios.  With  a  given  condition  number,  no  data  set  is  possible  for  which  a 
sidelobe  canceller  can  give  an  S/N  improvement  that  lies  above  this  curve.  Some  scenarios  marked 
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as  crosses  on  the  plot  are  illustrated.  The  most  interesting  data  sets  are  points  just  a  little  below 
the  curve,  representing  strong  jamming  that  can  be  deeply  nulled. 
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Figure  6-1.  Scenarios  and  performance  for  ideal  nulling  and  for  simulated  hardware 
parameter  choices. 


Individual  experiments  with  fixed  word-length  parameters  produce  reduced  S/N  improvements. 
Those  weights  obtained  by  accurately  solving  for  weights  using  an  inaccurately  computed  Cholesky 
matrix  are  marked  as  filled  circles  on  the  plot.  The  S/N  improvements  for  the  same  inaccurately 
computed  Cholesky  matrix  when  the  weights  are  computed  using  the  technique  of  section  8  are 
indicated  by  open  circles.  When  several  scenarios  with  the  same  theoretical  performance  were  tried, 
there  is  one  cross  and  a  scattering  of  filled  and  open  circles. 

The  data  in  Figure  6-1  are  for  the  final  parameter  choices. 

Suppose  we  are  experimenting  with  the  number  of  CORDIC  stages.  By  making  generous  choices 
for  the  other  finite  word-length  parameters,  wre  can  get  some  idea  of  how  parsimonious  we  could  be 
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with  the  number  of  CORDIC  stages.  The  other  finite  word-length  choices  were  guessed  at  similarly. 
Then,  by  varying  word-lengths,  guard  bits,  and  stages,  one  bit  or  stage  at  a  time,  we  could  zero  in 
on  the  best  trade-off  of  hardware  cost  and  performance. 

Several  observations  were  clear  from  studying  many  cases  in  this  way.  First,  the  loss  in  S/N 
improvement  due  to  parsimonious  parameter  choices  increases  with  the  condition  number  of  the 
data.  For  a  data  set  with  a  condition  number  around  300,  the  loss  is  negligible.  For  a  data  set  with 
condition  number  about  1000,  the  loss  is  about  2  dB.  Therefore,  if  we  want  a  system  to  be  capable 
of  50  dB  S/N  improvement,  it  would  be  unwise  to  try  to  make  parameter  choices  based  on  data 
sets  for  which  the  condition  number  is  higher  than  850.  Second,  when  one  parameter  is  chosen 
too  parsimoniously,  curves  that  plot  performance  as  a  function  of  another  parameter  can  be  too 
irregular  to  be  useful,  often  seeming  to  suggest  that  saving  a  bit  or  a  stage  is  helpful  rather  than 
harmful.  Third,  the  loss  in  S/N  improvement  does  not  seem  to  depend  on  the  number  of  jammers 
as  long  as  there  are  fewer  jammers  than  degrees  of  freedom.  All  the  experiments  illustrated  in 
Figure  6-1  were  run  using  35  CW  jammers. 

After  many  experiments,  final  parameter  choices  were  made.  The  number  of  CORDIC  stages, 
exclusive  of  the  three  stages  used  for  the  combination  of  CORDIC  correction  gain  and  forgetting 
factor,  was  selected  to  be  13,  e.g.,  vmax  —  12.  The  word-length  of  internal  registers  of  the  CORDIC 
cell  was  chosen  to  be  24  bits.  Of  these  24  bits  the  two  least  significant  bits  are  guard-bits  and 
are  not  passed  beyond  the  CORDIC.  The  tacked-on  vector  data  passed  between  CORDIC  cells  uses 
a  22-bit  word-length.  The  stored  Uj  are  also  22-bit  words.  Note  that  the  inputs  to  the  wafer, 
obtained  from  the  antenna  elements,  although  presented  as  22-bit  words,  must  not  be  allowed  to 
use  the  full  22-bit  dynamic  range  since  the  Cholesky  matrix,  generated  from  many  such  tacked-on 
vectors,  contains  the  energy  of  all  of  them  and  so  is  much  larger  than  any  single  one. 

It  was  decided  earlier  that  MUSE  should  be  capable  of  updating  its  Cholesky  matrix  about 
45,000  times  per  second.  Each  vector  needs  67  microcycles1  to  be  fed  into  MUSE,  so  the  time  for 
a  microcycle  must  be  about  1/3  microsecond.  The  system  clock  rate  is  four  times  as  fast,  about 
12  MHz. 

Statistical  accuracy  of  the  weights  to  within  about  1  dB  of  optimum  requires  about  320  up¬ 
dates  [3].  This  means  that  the  weights  can  be  updated  as  often  as  140  times  per  second.  However, 
MUSE  is  physically  capable  of  updating  its  weights  much  more  frequently  or  less  frequently  because 
a  “snapshot  cycle”  is  designated  by  a  special  external  signal. 


1  MUSE  accepts  a  new  64-element  vector  every  67  microcycles,  rather  than  every  65.  This  will 
be  explained  in  section  8. 
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7.  APPLICABILITY  OF  RESTRUCTURABLE  VLSI 


The  systolic  array  described  in  the  preceding  sections  has  32  supercells,  each  requiring  three 
CORDIC  cells.  A  CORDIC  cell  can  be  designed  as  an  integrated  circuit.  Such  a  circuit  requires  some 
tens  of  thousands  of  transistors,  classifying  it  as  VLSI.  Because  the  three  types  of  CORDIC  cells  are 
so  similar  to  one  another,  it  suffices  to  design  only  one  type  of  CORDIC  cell,  capable  of  operating  in 
three  different  modes.  The  MUSE  system  can  be  realized  efficiently  using  an  advanced  technology 
called  Restructurable  Very  Large  Scale  Integration  (RVLSI).  Lincoln  Laboratory  has  already  used 
this  technology  to  realize  six  different  systems,  so  the  technology  is  reasonably  mature.  The  RVLSI 
technology  is  briefly  described. 


7.1  Restructurable  Very  Large  Scale  Integration 


Restructurable  VLSI  comprises  a  design  methodology,  a  laser-based  interconnect  modification 
technology,  and  a  set  of  CAD  tools  for  building  large  area  integrated  circuits  [8,9].  Wafers  are 
fabricated  with  redundant  circuits  and  interconnect,  which  are  tested  after  fabrication.  A  laser 
is  then  used  both  to  form  and  to  break  connections  to  the  operable  circuits,  and  to  build  the 
desired  system.  The  laser  can  also  customize  circuitry  by,  for  instance,  setting  constant  coefficients 
in  circuits  used  to  implement  a  filter  function.  The  size  of  the  basic  replaceable  unit,  the  cell,  is 
determined  by  partitioning  considerations  and  fabrication  yield;  typically  a  cell  comprises  thousands 
of  transistors.  Experience  has  shown  that  building  a  wafer  with  about  twice  as  many  cells  as 
ultimately  needed  strikes  a  good  balance  between  interconnect  overhead  and  cell  yield. 

Several  laser  restructuring  technologies  have  been  developed  [10];  the  technique  to  be  used  in 
the  MUSE  application,  which  uses  the  laser  to  form  a  connection  between  two  adjacent  diffusions, 
is  completely  compatible  with  standard  IC  processing.  The  laser  diffused  link  is  used  for  both  signal 
and  power  connections.  Connections  are  broken  by  using  the  laser  to  vaporize  a  metallization  path. 
The  connections  and  the  metal  cuts  are  made  with  high  yield  and  appear  to  be  very  reliable. 

A  set  of  software  tools  has  been  developed  for  the  design  tasks  that  are  unique  to  wafer  scale 
design  [11].  Floorplanner  is  used  to  plan  the  tiling  of  a  wafer  with  cells  and  interconnect  and  IRT 
performs  a  unique  assignment  and  routing  based  on  test  results  for  each  wafer.  IRT  also  adds 
connections  for  testing  of  the  wafer  during  restructuring  and  allows  removal  and  replacement  of 
circuits  that  are  found  to  be  defective  during  this  process. 

Six  RVLSI  systems  have  been  built  on  three  different  wafer  designs  [8]  with  the  largest  one 
having  405,000  functional  transistors.  One  Integrator  wafer  scale  system  has  been  operating  in  a 
bench  tester  for  4  1/2  years  without  failure.  The  MUSE  system  and  another  now  being  built  [12] 
are  by  far  the  largest  and  most  logically  complex  wafer  scale  circuits. 
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7.2  Adapting  MUSE  for  RVLSI 


Several  design  choices  were  made  in  order  to  make  the  MUSE  design  even  more  compatible 
with  the  RVLSI  technology. 

First,  because  of  the  significant  cost  of  intercell  connections,  the  microcycle  was  divided  into 
four  steps  so  that  data  could  be  moved  into  and  out  of  a  CORDIC  cell  in  four  small  pieces.  Consider 
a  0  cell  that  must  read  two  22-bit  words  from  either  of  two  sources  (the  supercell  ahead  of  it  and 
the  supercell  behind  it)  while  outputting  two  22-bit  words  to  the  <f)  cells  of  its  supercell.  This 
would  have  required  132  connections  to  at  least  as  many  metallization  lines,  without  even  counting 
clocks,  control,  power  and  ground,  and  jumpers.  The  four-step  microcycle  allows  the  same  data  to 
be  moved  with  33  connections  and  33  metallization  lines. 

Second,  the  two  adder /subtractors  in  each  CORDIC  stage  were  replaced  by  one  adder /subtractor 
time-shared  between  two  tasks.  It  first  computes  its  “y”  output 

y  <-  y  +  ?pix, 

but  preserves  the  old  y  for  the  next  computation 
x+-x-  2lpiy. 

The  adder/subtractor  is  the  largest  part  of  the  stage,  much  larger  than  the  control  or  the  pipeline 
latches.  Therefore,  this  change  reduces  the  size  of  the  stage.  In  exchange,  it  takes  twice  a s  long  to 
perform  a  minirotation  since  this  is  essentially  the  time  to  perform  the  addition  and  the  subtraction. 
Since  a  CORDIC  cell,  even  after  this  change,  is  relatively  large,  and  since  an  adder  can  be  designed 
adequately  fast,  this  modification  was  judged  worthwhile. 

Third,  it  was  decided  to  provide  each  CORDIC  cell  with  adequate  memory  to  store  the  required 
real  or  imaginary  part  of  the  two  columns  of  the  Cholesky  factor,  even  though  the  0-CORDIC  cells 
have  no  need  for  this  memory.  The  memory  uses  only  a  minor  portion  of  the  silicon  area  of  each 
cell,  and  the  cost  of  this  silicon  area  must  be  balanced  against  the  competing  cost  of  providing 
connections  between  the  (fhce  11s  and  a  second  type  of  cell,  “memory,”  or  alternatively,  against  the 
cost  of  making  two  types  of  CORDIC  cells,  some  with  memory  and  some  without. 

The  fourth  change  is  very  much  more  complex  and  is  dealt  with  in  section  8.  Briefly,  after 
enough  samples  have  been  used  to  update  a  Cholesky  matrix  and  the  time  has  come  to  use  that 
matrix  to  solve  the  linear  equations  for  the  nulling  weights,  it  would  have  been  necessary  to  provide 
data  paths  to  pass  the  N(N  + 1)/2  matrix  elements  off  the  wafer  to  where  they  could  be  used  next. 
Instead,  we  found  a  way  to  use  the  CORDIC  cells  to  begin  the  process  of  solving  the  linear  equations. 
This  method  produces  2N  intermediate  results,  which  are  the  only  quantities  that  must  be  sent 
off  the  wafer,  and  they  are  moved  along  on  the  paths  already  provided  for  moving  the  tacked-on 
vector. 

Finally,  we  will  make  the  CORDIC  cells  in  two  mirror  image  versions  and  run  the  discretionary 
metallization  paths  between  them.  The  ability  to  use  CORDIC  cells  on  either  side  of  a  “bundle”  of 
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discretionary  connection  tracks  gives  extra  flexibility.  Details  of  the  RVLSI  CORDIC  cell  are  given 
in  section  9,  and  details  of  the  wafer  floorplanning  and  restructuring  are  given  in  section  10. 


7.3  Relationship  of  Cell  Yield  to  System  Yield 

In  this  subsection  it  is  assumed  that  the  process  by  which  a  wafer  is  manufactured  has  a 
significant  probability  of  defects.  Assume  we  make  K  cells  per  wafer.  Assume  that  each  cell  on  the 
wafer  functions  correctly  as  a  CORDIC  with  a  probability  pc,  independent  of  whether  its  memory 
works,  while  each  cell’s  memory  is  perfect  with  probability  pm ,  independent  of  whether  it  is  a  part 
of  a  functional  CORDIC  cell  or  not.  Therefore  a  cell  chosen  at  random  is  usable  as  a  0-CORDIC  with 
probability  pc  and  is  usable  as  a  ^CORDIC  with  probability  pcPm- 

Assuming  that  we  have  adequate  testability  and  connectability,  we  can  compute  the  probability, 
Pw ,  of  forming  a  completely  functional  MUSE  system  on  a  wafer  in  terms  of  AT,  pc,  and  pm.  A 
given  wafer  is  suitable  if  at  least  96  CORDIC  cells  work  and  if,  among  the  CORDlCs  that  work,  at 
least  64  have  working  memory.  The  probability  that  exactly  i  of  K  cells  have  CORDIC  function  is 


and  the  probability  that  exactly  i  working  CORDlCs  have  exactly  k  working  memories  is 


Pw  can  therefore  be  written  as 


This  formula  can  be  evaluated  by  a  computer,  but  it  is  useful  to  first  make  a  well-studied 
approximation  [7].  The  expected  number  of  working  CORDlCs  per  wafer  will  be  Nc  =  Kpc  and  its 
variance  will  be  =  Kpc(  1  —  pc).  Then 


The  approximation  is  valid  when  Kpc(  1  —  pc)  S>  1.  Both  the  approximation  and  the  correct 


formula  give  nonnegligible  results  except  when  i  is  within  about  five  standard  deviations  of  its 
mean. 


In  a  similar  fashion  the  expected  number  of  working  memories  among  i  CORDlCs  is  Nm(i)  = 
ipmi  with  variance  ipm(  1  —  pm).  Then  for  Pi/K.k/i  the  appropriate  approximation  is 
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(7.1) 


Pi/K-.k/i 


_  ( i-Kpc )2  _ (*  — »Pm) 

eXP  2Kpc(l-pc)  e  2'Pm(»-Pm) 
>/27r.firpc(l  -pc)  s/2nip^(T^Pm) 


This  approximation  is  valid  when  Kpc(  1  —  pc)  3>  1  and  kpm(  1  —  pm)  S>  1. 
Then  we  may  conveniently  compute  Pw  as 


«.=  E 


/t 

eXP  2/Cpcd-pc) 

V2lrKPc(l  -pc) 


t 


E 


_ ^Pm)2 

g  2tpm  ( 1  —  Pm ) 

v/27ripm(l  -  Pm) 


(7.2) 


A  computer  program  to  evaluate  this  double  sum  may  be  further  optimized  by  limiting  the 
summation  indices.  For  the  outer  sum  the  low  index  is  the  greater  of  96  and  Kpc  -  5y/Kpc(l  —  pc), 
and  the  upper  index  is  the  lesser  of  K  and  Kpc  +  5^Kpc(l  —  pc).  Similarly,  for  the  inner  sum,  the 
low  index  is  the  greater  of  64  and  ipm  —  5\/ipm(l  —  pm)  while  the  upper  index  is  the  lesser  of  i  and 
tPm  +  tyipm{l  “Pm)- 

In  Figure  7-1  the  solid  family  of  curves  give  wafer  yield,  for  a  wafer  design  using  130  cells, 
versus  CORDIC  cell  yield,  with  different  assumed  values  of  memory  yield.  The  curves,  from  right  to 
left,  are  for  memory  yields  of  50,  60,  70,  80  and  90  percent  (the  curves  for  80  and  90  percent  are 
too  nearly  identical  to  be  resolved). 

Of  course,  only  a  limited  number  of  wafers  will  be  made,  probably  ten.  It  is  not  difficult  to 
compute  Pfru,,  the  probability  that  the  best  of  the  ten  wafers  will  have  enough  useful  cells.  The 
probability  that  a  given  wafer  is  inadequate  is  1  —  Pw  and  the  probability  that  all  ten  are  inadequate 
is  (1  —  Pw )10.  Hence  the  probability  that  at  least  one  out  of  the  ten  wafers  can  be  used  for  a  wafer 
scale  MUSE  is 


Pine  =  1  -  (1  -  Pw)10- 

The  dotted  curves  in  Figure  7-1  give  the  probability  of  at  least  one  usable  wafer  versus  CORDIC 
cell  yield  for  the  same  choices  of  memory  yield  (curves  for  70,  80,  and  90  percent  are  almost 
indistinguishable) . 

An  examination  of  the  dotted  curves  in  Figure  7-1  shows  that  memory  failures  of  up  to  30  per¬ 
cent  make  almost  no  difference  in  the  expected  wafer  yield  or  probability  of  success.  Based  on  a 
CORDIC  yield  of  about  70  percent,  about  what  we  have  with  test  cells,  we  would  expect  about  an 
85  percent  chance  of  project  success  in  a  run  of  ten  wafers.  However,  this  probability  drops  sharply 
with  decreasing  CORDIC  yield. 

If  we  have  insufficient  CORDIC  cell  yield  to  make  a  complete  MUSE  system  on  a  single  wafer, 
we  can  still  make  a  wafer  scale  system  by  linking  together  two  wafers.  Only  a  few  connections 
need  to  be  made  between  the  two  wafers  because  a  relatively  small  number  of  connections  run 
from  supercell  to  supercell.  Roughly  speaking,  two  wafers  would  have  twice  as  many  CORDIC  cells 
to  start  with  and  therefore  should  allow  system  success  with  about  half  the  yield  that  would  be 
required  with  a  single  wafer.  In  reality  the  situation  is  slightly  more  favorable  since  the  lower 
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Figure  7-2.  Expected  system  yield  versus  CORDIC  cell  and  memory  yield,  using  best 
two  of  ten  wafers. 
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CORDIC  yield  implies  a  greater  variation  from  wafer  to  wafer  and  since  we  would  use  the  best  two 
wafers.  The  probability  of  forming  a  MUSE  system  from  the  best  two  wafers  out  of  ten  is  given  by: 


K  K  K  K 

P2w  =  EEE-E  S(n{ni,  n2,n3,...,ng,nio)) 

Til  =0  712=0  Tl3=0  7ljO=0 

p?(  1  -  pc)K~n >  ^  p?(  1  -  pc)K~n'  . . . 

(jr  \  n(ni,n-2 . nio)  /  \ 

p"10(i  -Pc)K_Bi°  E  r  )p™(i-Pmy-m 

”10/  m=64  W 


(7.3) 


where 


p(ni,n2,n3, . . .  ,ng,nio))  =  max(n,  +  n,), 


and 


S(m)  =  1  m  2  96 

[  0  otherwise 

In  Figure  7-2  we  plot  the  probability  of  finding  96  working  cordic  cells  with  64  working 
memories  on  the  best  two  out  of  ten  wafers.  The  curves,  from  right  to  left,  are  for  memory  yields  of 
50,  60,  70,  80,  and  90  percent  (curves  for  80  and  90  percent  are  indistinguishable).  From  Figure  7-2, 
we  would  predict  that  with  memory  yield  of  only  70  percent  and  CORDIC  yield  of  only  34  percent 
we  have  an  85  percent  chance  of  forming  a  MUSE  system  with  two  wafers  out  of  ten. 
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8.  USING  CORDIC  CELLS  TO  SOLVE  LINEAR  EQUATIONS  FOR 

WEIGHTS 


8.1  Mathematical  Preliminaries 


Consider  the  general  N  x  N  set  of  linear  equations 
AX  =  B, 


where  A  and  B  are  given  and  X_  is  to  be  determined.  To  solve  the  linear  equations,  we  will  introduce 
a  seemingly  unrelated  second  problem.  First,  construct  an  N  x  (N  +  1)  matrix  by  appending  B 
onto  the  right  edge  of  A. 


mi  a- 


Next,  postmultiply  it  by  an  (N  +  1)  x  (N  +  1)  unitary  matrix  Q  that  causes  the  last  column  of  this 
matrix  to  become  zeroes. 


[A\B}Q  = 


Third,  partition  Q  as 


Quu  I  Qt 


From  this  partition  we  can  obtain  the  equation 


AQr  +  qB  —  0 

leading  to 

-Q 

This  shows  that  the  solution  to  the  original  problem,  solving  for  A,  is  hidden  in  the  elements 
making  up  the  last  column  of  the  unitary  matrix  Q  in  the  second  problem.  If  we  can  zero  out  a 
column  B  tacked  on  to  a  matrix  A  using  a  unitary  transformation  Q,1  then  we  can  solve  AX  =  B. 
But  the  collection  of  CORDIC  cells  that  update  a  Cholesky  matrix  are  performing  just  such  a  unitary 

1  Actually,  there  is  no  requirement  that  Q  be  unitary. 
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transformation.  This  suggests  that  we  might  use  the  same  CORDIC  cells  to  solve  for  weights  W 
using  that  Cholesky  matrix  as  the  matrix  of  coefficients  in  the  linear  equations. 

Note  that  in  section  2  it  was  noted  that  there  are  generally  two  sets  of  linear  equations  to 
solve,  the  first  set  involving  L  and  the  second  set  involving  Lh.  But  in  the  case  of  a  sidelobe 
canceller  (a  particular  choice  of  the  steering  vector,  5  =  [0,0, ...  ,0, 1]*),  the  first  set  can  be  solved 
by  inspection  rather  than  by  computation,  so  that  only  the  second  set,  that  involving  Lh ,  needs  to 
be  solved  using  the  CORDIC  cells. 

The  first  set  of  equations  to  be  solved  was: 

hiVi  =  5i=0 

hiVi  +  I22IJ2  =  52  =  0 

him  +  h2V2  +  kzyz  =  53  =  0 


IniVi  +  •  •  •  +  InnVn  =  5  tv  =  1 

We  can  successively  see  that  3/1, 2/2*  -  -  -  *  VN- 1  are  0.  The  last  equation  in  the  set  reduces  to 
In,nVn  =  5tv  =  1, 

so  that  ?/TV  =  1/Inn-  In  short,  the  vector  Y  that  solves  LY  =  S  for  the  sidelobe  canceller  case  is 
Y  =  No  actual  computing  steps  are  needed — the  answer  is  available  by  inspection. 

However,  we  are  still  required  to  carry  out  computations  to  solve  for  ]V_  in  the  second  set  of 
equations  LhW_  =  H  =  j~S_-  Note  that  we  can  replace  the  scale  factor  by  any  other  which 
suits  us,  since  changing  it  would  only  scale  IT  by  a  constant,  which  cannot  affect  the  resulting 
signal-to-interference  ratio. 

In  the  following  discussion,  we  will  assume  that  in  each  supercell  the  two  columns  of  the 
snapshot  of  L  for  which  we  want  to  solve  for  \V  are  stored  in  a  snapshot  memory.  The  process 
of  solving  for  the  weights  can  therefore  be  interleaved  with  updating  the  Cholesky  matrix  as  more 
antenna  data  is  fed  into  MUSE. 

To  use  CORDIC  cells  to  solve  LhXV_  =  5,  we  begin  by  transforming  the  equation  set  into  one 
involving  a  lower- triangular  matrix  with  a  tacked-on  vector.  Define  the  “reversal”  matrix  J : 


J  = 


1 


1 


1 


1 
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If  we  premultiply  a  vector  or  matrix  by  J,  then  the  result  is  to  reverse  it  top  to  bottom.  If  we 
postmultiply  a  matrix  by  J,  then  it  is  reversed  left  to  right.  Note  that  JJ  is  an  identity  matrix. 
Now  we  manipulate  LhW_  =  5,  first  by  inserting  JJ  between  Lh  and  IT;  next  by  premultiplying 
by  J  to  give 


(JLhJ)(JW)  =  JS. 

And  finally  by  conjugating  the  entire  equation 
{JL'J^JWy  =  JS. 

This  equation  is  of  the  form  AX  =  B  where 
A  =  JLlJ , 

X  =  {JW)\ 

B  =  JS, 

and  the  form  of  A  is  a  lower  triangular  matrix: 


In,n 

In,n- i 

In-i, n-i 

In, 2 

In-i,2 

■  ■  h,2 

In,  i 

In-i,i 

■  ■  h,i  h,i 

8.2  Description  of  the  Q-operation 

Of  course,  if  we  solve  for  X  we  immediately  get  IT  since  IT  =  JX*.  The  tacked-on  vector 
B  is  [1,0,0,...  ,0]*.  Changing  Lh W  =  5  into  AX  =  B,  an  equation  involving  a  lower  triangular 
matrix,  has  been  purely  an  exercise  in  notation.  No  actual  computations  are  involved.  We  shall 
now  show  how  to  solve  AX  =  B  in  the  array  of  supercells  without  moving  the  Cholesky  matrix  out 
of  the  supercells  in  which  it  is  already  stored.  This  process  is  called  the  Q-operation. 

It  is  easy  enough  to  begin  the  process.  The  first  step  is  to  postmultiply  the  two  columns  made 
up  of  B  and  the  last  row  of  L  by  a  unitary  matrix  which  zeroes  out  the  leading  element  of  B. 

In,N  1 
In,n- i  0 

In, 2  0 

I.  In, i  0  . 
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lN^N  is  available  in  supercell  0,  where  it  was  created  in  its  guise  as  supercell  N-l.  That 
supercell  should  also  have  available  the  fixed  constant  &i  =  1.  The  supercell  determines  its  Co¬ 
operation  rotation  controls  (Q-op  controls)  in  both  its  0-CORDIC  and  its  (^CORDlCs  using  these 
two  quantities  as  column  leaders.  These  Q-op  controls  determine  the  operations  performed  on  the 
remainder  of  the  two  columns.  But  now  we  see  what  looks  like  a  problem.  The  next  step  involves 
the  pair  (In,n- i ,  0),  but  In,n- l  is  not  available  in  supercell  0.  It  is,  however,  available  in  supercell  1, 
which  already  has  a  connection  to  supercell  0.  We  can  pass  the  Q-op  controls  pi  from  supercell  0  to 
supercell  1  using  only  local  connections  and  then  carry  out  the  postmultiplication  which  originated 
with  the  pair  (In,N,  1)  in  supercell  0  on  the  corresponding  pair  {In,n- i,  0)  in  supercell  1.  Then  we 
pass  the  same  Q-op  controls  to  supercell  2  and  apply  them  to  the  pair  {In,N- 2,0).  We  can  move 
the  Q-op  controls  down  the  chain  of  supercells  into  supercell  31  in  this  way  and  then  back  up  the 
chain  all  the  way  back  to  supercell  0.  After  they  have  been  used  in  supercell  0  on  the  pair  (/#, i,  0), 
the  Q-op  controls  are  sent  off  the  wafer,  where  they  are  stored  for  later  use. 

As  a  result  of  all  this,  the  first  column  of  A  and  the  tacked-on  vector  B  are  modified.  No  further 
use  of  the  first  column  of  A  will  be  made.  Indeed,  we  don’t  need  to  save  it  as  it  is  computed.  The 
tacked-on  column,  f?,  has  been  modified  so  that  its  first  value  is  zero.  The  remaining  elements  of  i? 
(which  were  zero)  are  now  nonzero  values  [62,  63, . . . ,  b^Y  which  are  available  in  supercells  1,2,.. ., 
respectively,  where  they  were  created  by  the  rotations  just  discussed. 

The  next  step  is  to  postmultiply  the  second  column  of  A  and  the  modified  tacked-on  vector 
by  a  unitary  matrix  to  zero  out  the  next  element  of  the  modified  tacked-on  vector. 

In-x,n-i  fa 
In-1,N-2  fa 

In-  1,2  bx- 1 

„  In-  1,1  bx  m 

Both  elements  making  up  the  first  pair,  62),  are  available  in  supercell  1,  62  having 

just  been  computed  there.  These  are  the  leaders  used  to  set  up  Q-operation  rotation  controls  in 
supercell  1.  These  Q-op  controls  are  moved  to  supercell  2  where  they  control  the  CORDlCs  acting 
on  the  pair  (/Ar-i,Ar_2, 63),  then  on  to  supercell  3  where  they  control  the  CORDlCs  acting  on  the 
pair  (In-\,N-3,  fa),  etc.  Ultimately  these  Q-op  controls  also  turn  around  at  supercell  31  and  travel 
up  the  chain  and  are  sent  off  the  wafer  at  supercell  0.  These  Q-op  controls  are  also  stored  for  later 
use.  At  the  end  of  this  step  the  elements  of  the  now  twice  modified  tacked-on  vector  [63,64,...]* 
are  available  in  supercells  2,3, . . .,  exactly  where  they  will  be  needed  to  carry  out  the  third  step. 
(The  modified  second  column  of  A  is  no  longer  needed.) 

A  step  during  which  Q-op  controls  are  determined  by  vectoring  is  called  a  master  Q-operation . 
The  steps  during  which  these  angles  are  used  in  CORDIC  cells  downstream  are  called  slave  Q - 
operations. 

Using  this  procedure,  we  can  carry  out  the  entire  process  of  zeroing  out  the  tacked-on  vector 
using  CORDIC  rotations  of  data  that  is  never  moved  from  supercell  to  supercell.  The  only  information 
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that  moves  during  the  Q-operation  is  the  Q-op  rotation  controls.  These  are  passed  only  locally,  to 
the  adjacent  supercells,  so  no  global  communication  paths  are  needed.  Each  set  of  Q-op  controls  is 
ultimately  passed  off  the  wafer,  where  they  are  to  be  used  in  the  final  phase  of  weight  determination. 

In  this  final  phase,  the  elements  in  the  last  column  of  the  unitary  matrix  Q  must  be  determined. 
Although  we  have  used  the  cordics  to  postmultiply  [A  |  B]  by  Q,  Q  was  never  determined  in  the 
standard  matrix  form.  The  Q-op  controls  which  have  been  passed  off  the  wafer  are  all  that  is 
needed  to  instruct  CORDIC  cells  to  carry  out  the  operations  which  realize  a  postmultiplication  by 
Q .  Therefore  one  can  use  any  CORDIC  cells  with  these  same  controls  to  postmultiply  an  identity 
matrix  by  Q .  This  will  give  the  elements  of  Q,  which  can  be  used  directly  as  the  desired  weights  in 
the  nulling  problem. 


8.3  Summary  of  Revised  MUSE  Control  and  Timing 

The  systolic  array  presented  in  earlier  sections  was  100  percent  efficient.  Therefore  there  can 
be  no  time  available  for  the  CORDlCs  to  carry  out  the  Q-operation.  It  was  decided  to  interleave  the 
steps  of  the  Q-operation  with  subsequent  updating  of  the  Cholesky  factor,  by  providing  two  “idle” 
microcycles  in  every  macrocycle.  It  is  also  necessary  to  revise  the  latency  to  compensate  for  the 
longer  macrocycle.  Here  the  revised  timing  is  summarized. 

The  fastest  clock  present  in  the  system  is  called  a  half-step .  A  half-step  is  the  time  needed 
to  pass  11  bits  between  CORDIC  cells.  Two  consecutive  half-step  intervals  make  a  step  and  suffice 
to  move  a  word.  The  reciprocal  of  the  step  interval  is  the  rate  at  which  new  data  is  presented  to 
an  adder  within  any  CORDIC  stage.  Two  steps  are  called  a  microcycle.  A  microcycle  is  the  time 
required  to  input  one  complex  element  into  any  0-CORDIC,  and  its  reciprocal  is  the  rate  at  which 
new  rotations  may  begin  in  any  CORDIC  cell. 

The  time  interval  between  input  of  successive  64-element  vectors  into  MUSE  is  67  microcycles, 
called  a  macrocycle .  At  the  input  to  any  0-CORDIC  cell  the  67  microcycles  are  assigned  to  move 
data  as  described  in  section  5  except  that  a  blank  microcycle  is  inserted  between  each  tacked-on 
vector.  These  blank  microcycles  are  used  to  move  Q-op  controls. 

Any  macrocycle  can  be  designated  as  a  “snapshot”  cycle.  When  the  Cholesky  factor  has  been 
updated  to  reflect  the  effect  of  the  vector  fed  into  MUSE  during  a  snapshot  cycle,  the  Q-operation 
is  initiated.  The  interval  between  successive  snapshot  cycles  is  called  an  update  cycle. 

The  timing  of  MUSE  all  builds  up  from  the  half-step.  The  CORDIC  cell  designed  is  capable  of 
being  clocked  at  12  MHz.  Table  8-1  summarizes  the  various  clock  speeds  and  periods  implied  by 
the  12  MHz  half-step  clock. 

Since  MUSE  has  96  CORDIC  cells,  each  performing  three  million  new  rotations  per  second,  a 
conventional  computer  requiring  ten  instructions  per  rotation  would  have  to  perform  2.88  billion 
instructions  per  second  to  perform  MUSE’s  task. 
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The  revised  latency  of  a  supercell,  in  units  of  microcycles,  must  be 


TABLE  8-1. 

Projected  Speed  of  MUSE 


Interval 

Duration 

Rate 

Half-Step 

83  1/3  ns 

12  MHz 

Step  =  2  half-steps 

166  2/3  ns 

6  MHz 

Microcycle  =  2  steps 

333  1/3  ns 

3  MHz 

Macrocycle  =  67  microcycles 

22  1/3  /is 

44.78  kHz 

assuming 

Update  Cycle  300 

macrocycles 

6.7  ms 

150  per  second 
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9.  RVLSI  CELL  DESCRIPTION 


The  CORDIC  cell  has  been  fabricated  in  2 fim  CMOS  and  tested  at  the  design  speed.  It  amounts 
to  54,000  transistors  and  is  5.5  mm2  x  5.6  mm2  excluding  output  drivers  and  bond  pads.  This 
section  briefly  describes  the  cell  layout.  It  consists  primarily  of  a  large  datapath,  two  static  RAMs, 
and  a  chunk  of  control  logic. 

9.1  Cell  Datapath 

The  main  datapaths  and  functional  blocks  of  the  CORDIC  cell  are  shown  in  Figure  9-1.  Eleven- 
bit  data  comes  from  the  forward  and  reverse  connections  at  the  top.  It  is  turned  into  22-bit  words 
and  passed  into  the  input  switching  section,  then  down  through  the  CORDIC  stages  and  output 
switching  before  being  split  into  11-bit  words  again.  The  delay  memory  is  provided  for  the  <f>- 
CORDics  to  delay  L{\  the  snapshot  memory  stores  one  complete  L  for  use  in  the  Q-operation. 
Along  the  right  edge,  the  intracell  serial  communication  path  is  shown;  this  carries  information 
during  the  Q-operation. 

Inside  the  input  switching  box  are  several  delay  registers,  two  shift  registers,  two  adders,  and 
a  big  multiplexer.  The  shift  registers  receive  serial  data  from  the  serial  path  and  can  insert  it 
into  the  main  datapath;  this  function  is  used  by  the  0-CORDIC  cells  for  the  Q-operation.  The  gain 
compensation  is  split  into  three  adders;  the  first  two  are  in  the  input  section  to  scale  down  the  input 
number  slightly  (by  the  factor  138/256).  This  ensures  that  preventable  overflows  are  avoided. The 
multiplexer  allows  one  of  six  of  data  values  to  be  entered  into  the  datastream  going  to  the  CORDlCs: 
the  input,  the  memory  output,  either  serial  register,  the  constant  0,  or  the  constant  218.  The  last 
two  are  only  used  by  the  Q-operation.1 

The  basic  structure  of  the  CORDIC  stage  was  described  in  section  3;  the  .sharing  of  one  adder 
between  real  and  imaginary  parts  of  a  word  is  described  in  section  7.2.  The  resulting  logic  consists 
of  three  registers,  a  one-of-two  multiplexer,  and  an  add/subtract  unit.  The  CORDIC  stages  are  24 
bits  wide;  13  stages  are  required  to  meet  the  accuracy  criterion  described  in  section  6. 

Within  the  output  switching  box  are  several  delay  registers,  one  serial  register,  and  one  adder. 
The  serial  register  puts  data  onto  the  internal  serial  path  during  the  Q-operation  in  the  </>-  CORDlCs. 
The  adder  completes  the  gain  compensation  (with  a  gain  of  9/8).  The  delay  memory  simulates  a 
long  shift  register  with  a  circular  buffer.  Address  generation  is  provided  by  a  simple  counter  that 
resets  to  0  upon  reaching  51;  each  location  is  first  read,  then  written.  The  snapshot  memory  stores 
two  L  vectors  acquired  at  different  times  as  outlined  in  section  5.  Its  address  generation  is  more 
complicated  and  is  described  in  the  next  subsection. 


1  Since  the  input  and  output  word-lengths  of  a  CORDIC  cell  are  the  same,  overflow  cannot  be 
made  impossible.  For  example,  if  the  real  and  imaginary  inputs  are  approximately  equal  and  the 
angle  of  rotation  is  near  to  -45°,  we  expect  the  real  part  of  the  output  to  exceed  either  input  by 
about  \2.  While  overflow  for  all  possible  inputs  cannot  be  prevented,  we  can  insist  that  overflows 
occur  only  when  the  correct  output  data  would  not  be  representable  in  a  22-bit  word. 
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Figure  9-1.  Datapaths  of  CORDIC  computation  chip. 
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The  serial  path  is  used  for  a  number  of  purposes,  depending  on  whether  the  cell  is  a  6  cell 
or  a  0  cell.  In  all  cases  it  is  only  used  during  the  Q-operation.  Q-op  control  data  is  passed  from 
supercell  to  supercell  during  the  “blank,”  but  within  the  supercell  the  blank  contains  rotation 
data.  Therefore,  the  6  cell  receives  parallel  format  angle  values  and  passes  them  serially  to  the  <p 
cells.  The  (f>  cells  must  output  the  angle  value  in  parallel  format.  During  a  slave  Q-operation,  each 
CORDIC  stage  must  read  one  of  the  angle  bits,  while  during  a  master  Q-operation  each  stage  must 
set  one  of  them.  The  internal  serial  path  provides  all  of  this  capability;  the  control  section  ensures 
that  the  operations  occur  at  the  right  time. 


9.2  Control 

MUSE  is  a  systolic  array  and  therefore  depends  on  local  communication  for  control.  As  men¬ 
tioned  in  earlier  sections,  there  are  two  control  signals:  D  for  direction  control,  and  Q  to  indicate 
when  a  snapshot  should  occur.  Q  may  reach  a  6  cell  from  either  the  forward  or  the  reverse  direc¬ 
tion.  All  control  is  generated  from  these  two  inputs,  plus  one  global  input,  ZERO.  One  update  cycle 
consists  of  the  following:  a  snapshot  in  the  forward  direction,  a  snapshot  in  the  reverse  direction, 
N  —  2  slave  Q-operations,  and  two  master  Q-operations.  The  Q-operations  occur  in  both  forward 
and  reverse  directions  and  may  be  interleaved. 

The  outside  world  designates  the  beginning  of  a  forward  snapshot  by  asserting  the  Q  line  to 
the  first  cell  in  the  chain.  The  particular  column  of  L  next  updated,  called  a  snapshot  vector, 
will  be  operated  on  by  each  cell,  and  then  copied  into  the  snapshot  memory.  Each  cell  sees  this 
snapshot  vector  in  two  pieces,  first  in  the  forward  direction  and  then  in  the  reverse  direction.  Each 
part  of  the  vector  is  a  different  length  as  described  in  section  5,  but  the  sum  of  the  two  lengths  for 
a  given  cell  is  always  N+l.  The  forward  vector  is  copied  into  the  snapshot  memory  from  address 
0  upwards,  while  the  reverse  vector  is  copied  into  the  memory  from  address  N  downward. 

Each  supercell  in  the  chain  follows  the  same  rules  for  the  Q-operation:  when  a  forward  Q- 
operation  is  received,  use  the  data  pointed  to  by  the  reverse  pointer  and  decrement  the  pointer. 
If  the  reverse  pointer  is  0,  generate  a  forward  master  Q-operation  in  the  next  blank.  When  a 
reverse  Q-operation  is  received,  use  the  data  pointed  to  by  the  forward  pointer  and  increment  the 
pointer.  If  the  forward  pointer  is  N,  generate  a  reverse  master-Q  operation  in  the  next  reverse 
blank.  Whenever  a  master  Q-operation  is  generated,  pass  a  slave  Q-operation  on  down  the  chain. 
Always  pass  along  a  slave  Q-operation. 

Once  the  reverse  snapshot  hats  finished,  Q-operations  may  occur.  In  the  first  supercell,  the 
reverse  pointer  is  0  because  its  forward  snapshot  is  N  words  long.  Therefore,  it  generates  a  forward 
master-Q  operation  as  soon  as  the  reverse  snapshot  is  done,  and  passes  a  slave  Q-operation  to 
the  second  supercell.  The  second  supercell,  whose  reverse  pointer  was  1,  decrements  it  to  0  and 
then  generates  a  master  Q-operation.  This  pattern  continues  down  the  chain  and  back  up  again, 
so  that  the  first  supercell  finally  carries  out  63  slave  Q-operations  followed  by  its  second  master 
Q-operation.  One  Q-op  control  value  is  passed  off  the  wafer  in  the  reverse  blanks  for  each  reverse 
Q-operation  the  first  supercell  carries  out. 
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10.  WAFER  SCALE  CONNECTIONS 


MUSE  wafer  wiring  is  simply  an  iteration  of  the  connections  for  one  supercell  shown  in  Fig¬ 
ure  10-1,  plus  certain  peripheral  circuits  and  connections.  Note  that  there  are  three  general  classes 
of  connection  nets.  The  first  comprises  data  and  the  associated  D  and  Q  control  signals,  a  bundle 
of  11  to  13  conductors;  the  second,  the  serial  paths  between  rotator  cells  of  a  supercell;  and  the 
third,  global  clock  and  control  signals,  plus  power.  The  bundled  and  serial  classes  together  can  be 
described  as  “interconnections,”  in  contrast  to  global  connections.  The  RVLSI  CORDIC  cells  have 
interconnection  ports  along  one  edge,  global  ports  along  the  opposite  edge.  They  are  placed  as  mir¬ 
rored  images  between  alternating  types  of  wiring  channels,  corresponding  to  the  interconnections 
and  the  global  connections. 

10.1  Interconnection 


The  bundled  interconnections,  which  are  a  realization  of  the  local  communication  of  Figure  5-2, 
are  illustrated  in  a  simplified  form  in  Figure  10-2.  The  signal  paths  are  serpentine,  with  a  left-to- 
right  flow  in  one  channel,  right-to-left  in  the  next,  with  each  bundle  represented  as  a  single  track. 
The  figure  shows  the  interconnections  of  eight  supercells,  indicated  by  the  subscripts,  assembled 
from  24  CORDIC  cells.  As  on  a  real  wafer,  assignment  is  constrained  by  cell  yield.  Some  cells,  such 
as  those  marked  with  central  crosses  in  Figure  10-2,  cannot  be  used,  and  others,  with  defective 
memories  corner  crosses  can  be  used  only  for  the  6  function. 

Note  that  in  the  upper  channel  the  “F,”  “R,”  and  “OUT”  ports  are  in  the  order  of  forward 
signal  flow;  if  one  has  complete  freedom  of  cell  assignment,  as  for  supercells  0  and  1,  only  three 
tracks  are  necessary.  In  the  lower  channel  the  port  order  is  contrary  to  signal  flow,  and  four  tracks 
are  required,  a s  for  supercells  4  and  5.  In  order  to  use  cells  with  defective  memory  as  6  cells, 
it  is  often  necessary  to  use  an  extra  track,  as  shown  for  supercells  2  and  3  of  the  forward-flow 
channel  and  6  and  7  of  the  reverse-flow  channel,  so  that  as  many  as  five  tracks  may  be  required. 
In  experiments  with  a  number  of  randomly  placed  defective  cells,  the  IRT  routing  tool  [11]  never 
required  more  than  five  tracks,  so  this  can  be  taken  as  a  practical  upper  limit. 

Actual  interconnect  is,  of  course,  for  bundles,  not  single  tracks.  The  interconnect  channel  must 
accommodate  five  13-track  bundles,  as  many  as  five  more  tracks  for  serial  signals,  a  test  track  for 
proof  of  restructuring  connections,  and  spares  for  possible  defective  tracks.  Therefore,  a  channel 
of  85  tracks  was  chosen,  with  the  connections  to  the  CORDIC  cells  by  an  array  of  five  link  fields  of 
the  type  shown  in  Figure  10-3,  with  potential  links  shown  as  circles.  These  links  can  connect  the 
tracks  to  cell  ports,  above  and  below  the  channel,  via  short  vertical  tracks  known  as  stubs.  The 
stub  labeling  is  consistent  with  Figure  10-1  except  for  the  stubs  labeled  XIN ,  XIN2 ,  and  XOO , 
which  are  spare  ports  for  SIN ,  SIN 2,  and  SOO. 

The  top  track  (labeled  “16”)  is  for  connection  testing,  and  the  remaining  16  are  for  interconnect. 
The  data,  Q,  and  D  signals  have  links  in  a  pair  of  chevrons,  spaced  three  tracks  apart.  Primary 
links  for  these  13  signals  are  in  the  lower  chevron,  spares  in  the  upper.  The  final  three  spares  are 
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Figure  10-1.  MUSE  supercell  connections. 
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10.1.1  Channel  Defect  Capacity 


Each  signal  track  can  quite  naturally  be  assigned  to  one  of  three  classes,  according  to  the 
modulo  3  value  of  its  position  in  a  17-track  subarray,  as  in  Table  10-1. 

A  worst-case  connectivity  requirement  uses  five  bundles  of  tracks  per  channel.  Without  count¬ 
ing  the  tracks  used  for  serial  signals  within  a  supercell,  this  will  use  up  25  tracks  in  class  0,  out 
of  thirty  available,  and  20  tracks  in  each  of  classes  1  and  2,  with  25  of  each  available.  Therefore  a 
channel  may  have  up  to  five  defective  tracks  per  class.  However,  after  consideration  of  the  serial 
signal  track  requirements,  a  channel  must  not  have  more  than  a  total  of  ten  defective  tracks. 


10.1.2  Assignment  and  Routing 

Assignment  is  the  selection  of  particular  RVLSI  cells  to  serve  as  specified  processors;  routing 
is  the  selection  of  specific  nets  of  track  segments  to  interconnect  the  assigned  cells.  Assignment 
for  MUSE  is  straightforward,  with  cell  order  pretty  much  following  the  sequence  of  the  logical 
description  of  the  wafer.  However,  to  minimize  track  requirements  it  is  important  in  the  “contra¬ 
flow”  channels  to  have  the  0  cell  of  a  supercell  to  the  left  of  its  (j)  cells. 

A  MUSE  wafer  can  be  routed  with  a  minimum  number  of  tracks  if  a  specific  algorithm  is  used. 
The  basic  idea  is  to  proceed  in  the  sequence  the  tracks  appear  in  Figure  10-3,  bottom  to  top:  All 
data  bit  “0”  nets  are  routed  first,  then  data-bit  “1,”  and  so  forth  through  the  Q  and  D  nets.  For 
each  signal,  primary  tracks  are  used  first;  secondary  ones  are  used  only  if  a  primary  is  not  available. 
Also  for  each  signal  the  longest  nets  are  joined  first,  i.e.,  those  with  four  ports,  then  three,  and 
finally  two.  Serial-signal  nets  are  joined  last,  again  with  longer  ones  preceding  shorter.  If  no  track 
linking  to  an  “S”  stub  is  available  the  corresponding  “X”  stub  may  be  used  instead,  as  they  are 
connected  within  the  cell. 


10.1.3  Global  Connections 

Global  connections  are  much  simpler  than  interconnections  because  every  active  cell  is  con¬ 
nected  in  almost  exactly  the  same  way.  Restructuring  links  are  used  for  three  purposes.  First, 
the  power  connections  are  supplied  with  power-decoupling  capacitors,  which  are  separately  tested. 
Those  not  shorted  are  connected  to  the  Vdd  bus  by  “zapping”  large  (200  /zm,  4  ohm)  links.  Second, 
there  are  separate  power  supply  connections  for  the  cell  memories  and  the  data  path.  Third,  there 
are  spare  tracks  available  for  the  global  signals  ZERO ,  REAL ,  CLK1F ,  and  CLK2F .  1 


10.2  Testing  During  Restructuring 

Two  types  of  testing  are  concurrent  with  wafer  restructuring.  The  first  is  a  continuity  test  for 
link  formation.  Every  CORDIC  cell  port  has  an  associated  photodiode.  When  an  interconnection 

1  These  power  links  are  placed  within  the  cells. 
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TABLE  10-1. 


Class  Assignment  of  Tracks  for  Defect  Avoidance 


Class 

Bundled  Signals 

Tracks  per  Channel 

0 

Data  0,3,6, 9;  S 

30 

1 

Data  1,4,7,10 

25 

2 

Data  2,5,8;  S 

25 

net  is  formed,  one  of  the  stubs  of  the  net  is  also  linked  to  a  test  track,  number  16  in  the  link  field 
of  Figure  10-3,  which  is  in  turn  connected  to  a  wafer-edge  pin.  The  photodiodes  of  the  net  can 
then  be  consecutively  illuminated  by  a  low-power  laser  beam  while  photocurrent  is  sensed  at  the 
pin.  Upon  completion  of  the  test  the  stub  is  cut  away  from  the  test  track. 

The  second  test  is  of  functionality  of  each  cell  as  it  is  added  to  the  circuit.  For  the  MUSE 
wafer,  this  means  adding  to  the  circuit  the  three  CORDIC  cells  of  a  supercell,  and  then  applying 
stimuli  to  this  added  supercell  while  observing  signal  outputs.  The  technique  is  illustrated  in  Figure 
10-4,  which  shows  a  step-by-step  formation  of  the  eight  supercell  simplified  ‘‘wafers”  of  Figure  10-2. 

Step  A  Three  test  bundles  (shown  as  tracks)  are  linked  up  to  provide  access  to  all 
supercell  inputs. 

Step  B  The  last  supercell  of  the  chain  (supercell  7  in  this  example)  is  connected,  with 
its  output  and  reverse  input  brought  out  to  wafer-edge  pins,  and  its  forward 
input  connected  to  one  of  the  test  bundles.  This  single  supercell  is  tested. 

Step  C  The  next-to-last  supercell  is  connected  and  tested  in  conjunction  with  the  first, 
with  forward  input  connected  to  a  different  test  bundle.  The  test  bundle  used 
in  Step  B  is  now  re-used  as  a  permanent  connection  from  the  outputs  of  the 
two  </>-CORDlC  cells  of  supercell  6. 

This  procedure  is  repeated  for  each  supercell  in  turn,  until  all  have  been  connected.  The  important 
point  is  that  the  test  bundles  use  only  track  segments  that  will  eventually  be  used  in  the  interconnect 
itself.  At  each  stage  in  Figure  10-4,  new  interconnect  is  formed  by  cutting  off  segments  of  the  old 
test  bundle.  The  only  extra  wafer  wiring  resources  required  are  an  additional  set  of  wafer-edge  pins 
at  the  beginning  of  the  chain. 
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Figure  10-2.  Symbolic  bundled  interconnections  of  eight  MUSE  supercells  in  two 
wiring  channels . 
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Figure  10-3.  Link  field  for  an  interconnect  channel.  The  stubs  are  labeled  for  their 
eventual  cell  ports  above  or  below  the  channel. 
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Figure  10-4.  Step-by-step  connection  of  supercells  of  a  MUSE  wafer.  The  active  test 
track  at  each  step  is  marked  by  a  heavy  line;  inactive  ones  by  light  lines. 
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11.  MUSE  INVESTIGATION  TEST  EQUIPMENT  (MITE) 


11.1  Overview 


A  test  and  demonstration  system  (MUSE  Investigation  Test  Equipment)  has  been  developed 
to  permit  us  to  investigate  and  demonstrate  the  performance  of  the  MUSE  wafer.  MITE  can  send 
inputs  to  MUSE  at  full  speed  and  collect  the  angle  data  representative  of  the  resulting  weights. 
Extensive  software  permits  MITE  to  generate  test  data  sets,  interpret  and  display  the  results  as 
antenna  patterns  or  nulling  S/N  improvement,  and  compare  the  results  with  results  of  a  bit-by-bit 
simulation  of  the  wafer. 

The  primary  objective  of  MITE  is  to  test  the  full  wafer.  But  prior  to  actual  wafer  fabrication, 
MITE  has  been  used  to  test  individual  CORDIC  cells  in  a  subscale  processor  that  will  be  discussed 
in  appendix  B.  This  subscale  test  is  important  for  several  reasons: 


•  It  verifies  the  correct  operation  of  MITE  software  and  hardware. 

•  It  verifies  the  correct  operation  of  multiple  CORDIC  cells  operating  together  in 
the  MUSE  configuration  and  with  the  clocks  and  controls  of  the  final  wafer 
scale  system. 

•  It  verifies  the  concept  of  operating  a  subset  of  supercells  in  the  manner  in  which 
they  will  be  tested  during  the  restructuring  of  the  MUSE  wafer,  as  described 
in  section  10.3. 


The  subscale  tests  for  N  =  8,  using  12  CORDIC  cells,  have  verified  all  aspects  of  the  system 
designs  of  both  MUSE  and  MITE. 

MITE  allows  an  investigator  to  create  sets  of  N-element  data  vectors  drawn  from  a  common 
statistical  ensemble,  representing  signals  that  could  have  been  observed  given  a  specific  antenna 
geometry  and  jammer  scenario.  Once  the  data  sets  have  been  created,  they  are  stored  in  a  memory, 
from  which  they  may  be  sent  into  MUSE  at  a  high  speed  in  any  order  specified  by  the  investigator. 
The  investigator  may  create  several  different  data  sets,  representing  different  jammer  arrangements, 
to  examine  the  performance  of  MUSE  as  the  statistics  of  the  input  data  changes  with  time.  MITE 
is  also  able  to  initiate  snapshot  cycles  in  MUSE  and  obtain  the  resulting  Q-operation  angle  controls, 
which  are  captured  in  a  second  memory.  The  MITE  then  converts  the  Q-operation  angle  controls 
into  a  corresponding  weight  vector  and  can  generate  graphical  displays  of  the  resulting  performance 
of  MUSE. 

Figure  11-1  shows  the  MITE  test  system  with  a  four-supercell  subscale  MUSE  configuration 
under  test. 
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Figure  11-1.  Photograph  of  MUSE  Investigation  Test  Equipment  (MITE). 

11.2  MITE  Architecture 

Figure  11-2  shows  a  block  diagram  of  the  MITE.  A  PC/AT  computer  is  used  as  the  host  of  the 
system.  The  PC/AT  is  responsible  for  generating  input  data,  converting  Q-operation  angle  controls 
to  weight  vectors  and  generating  displays  and  performance  measures.  The  interface  system  consists 
of  two  large  memories  and  their  respective  controls.  The  interface  system  has  the  responsibilities  of 
supplying  MUSE  with  complex  data  samples,  initiating  snapshot  cycles,  and  gathering  the  resulting 
weight  vector  information  from  MUSE. 


11.3  Input  Data  Memory  and  Control 

The  data  input  system  has  three  control  signals  to  MUSE:  D ,  Q,  and  ZERO.  Their  respective 
functions  are  described  in  detail  in  section  9.  The  PC/AT  uses  an  interrupt  once  per  macrocycle 
from  the  interface  system  to  control  Q  and  ZERO.  The  direction  control,  D,  is  generated  by  the 
control  logic. 
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SYSTEM  HOST  SYSTEM  INTERFACE  MUSE  SYSTEM 


Figure  11-2.  MUSE  Investigation  Test  Equipment. 


The  input  memory  consists  of  a  32k  by  48-bit  static  memory.  This  memory  is  divided  into 
512  vectors  of  64  complex  elements.  Each  memory  location  corresponds  to  one  complex  element 
of  a  vector.  The  upper  nine  bits  of  the  address  specifies  an  entire  vector,  and  the  lower  six  bits 
correspond  to  a  complex  element  within  that  vector.  Both  the  PC/AT  host  and  the  MUSE  interface 
system  have  access  to  the  memory.  The  PC/AT  is  able  to  take  control  of  the  memory’s  data  and 
address  busses  and  write  data  vectors  into  the  appropriate  sections  of  the  memory. 

The  addressing  of  the  memory  by  the  interface  hardware  is  controlled  by  a  DSP56001A  mi¬ 
croprocessor  and  a  counter.  The  counter  generates  the  lower  six  bits  of  the  address,  which  must 
be  incremented  every  microcycle.  The  DSP56001A  is  used  to  generate  the  upper  nine  bits  of  the 
address,  which  may  be  changed  in  a  more  intricate  manner,  but  only  once  per  macrocycle.  During 
initialization,  information  is  sent  from  the  PC/AT  to  the  DSP56001A  that  specifies  the  location  of 
each  data  set  in  the  memory,  the  order  in  which  the  data  sets  are  to  be  used  and  the  number  of 
vectors  to  use  from  each  data  set. 


11.4  Output  Data  Memory  and  Control 


The  output  memory  and  buffers  are  shown  in  Figure  11-3.  The  latches  latch  the  four  11-bit 
output  values  from  MUSE  during  a  microcycle  in  which  output  data  is  present,  indicated  by  the  Q 
strobe  from  MUSE.  Then  at  a  much  slower  rate,  the  data  is  written  into  the  FIFO.  A  counter  is 
used  to  count  the  write  strobes  and  interrupt  the  PC/AT  when  a  complete  set  of  angle  words  has 
been  accumulated  and  is  ready  for  the  PC/AT  to  read. 

The  output  data  control  signals  are  generated  in  set  of  PALs. 


61 


144244-15 


Figure  11-3.  Output  memory  and  control. 


11.5  User  Interface  Software 


11.5.1  Pull-down  Menus 

The  software  system  is  completely  menu-driven.  The  system  allows  the  investigator  to  create, 
edit,  view,  and  load  the  different  types  of  parameter  files  used  by  the  system.  These  parameter  files 
are  stored  on  the  hard  drive  of  the  PC/AT  that  allows  the  parameter  files  to  be  used  at  a  latter 
date. 


11.5.2  Scenario  Generation 

The  investigator  may  create  jamming  scenarios  in  two  ways.  The  first  type  of  data  set  is 
based  on  antenna  patterns.  The  investigator  can  create  a  parameter  file  describing  a  simple  main- 
beam  antenna,  specifying  its  beam  width  and  sidelobe  levels,  and  another  parameter  file  with  the 
locations  and  gains  of  the  omnidirectional  auxiliary  antennas,1  and  yet  a  third  parameter  file  that 
specifies  the  directions  and  powers  of  multiple  jammers.  Using  these  three  parameter  files,  the 
PC/AT  can  create  a  scenario  file,  which  is  a  collection  of  vectors  drawn  from  a  population  with 
statistics  appropriate  to  the  jamming  and  antenna  geometry  described. 

Alternatively,  the  investigator  may  create  a  scenario  file  with  no  necessary  physical  significance, 
but  well  suited  to  testing  the  fundamental  S/N  improvement  limitations  of  MUSE  without  regard 

1  These  may  be  specified  individually  or  placed  randomly  by  the  PC/AT. 
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to  the  specifics  of  the  antenna  geometry  or  jammer  locations.  The  mathematical  description  of  this 
type  of  scenario  is  described  in  appendix  A. 

11.5.3  System  Operation 

Scenario  files  must  be  loaded  into  the  input  memory.  The  investigator  may  select  any  desired 
scenario  files  from  the  list  of  all  previously  created  scenarios  files,  using  the  menu  system.  Next 
the  investigator  loads  an  addressing  information  file,  containing  information  that  tells  the  system 
in  what  order  to  use  the  scenario  files,  and  how  many  vectors  to  use  from  each  set. 

The  investigator  must  then  load  a  file  containing  a  list  of  snapshot  cycles.  These  are  the 
macrocycles  in  which  the  system  will  order  MUSE  to  start  the  process  to  solve  for  the  weight 
vector.  The  system  will  produce  a  weight  vector  for  each  snapshot  cycle  in  the  file. 

The  investigator  is  then  able  to  start  the  system  and  then  view  the  results.  The  setup  infor¬ 
mation  and  resulting  weight  vectors  are  written  to  a  results  file. 

The  performance  of  the  weight  vectors  may  be  displayed  in  two  different  formats.  The  first 
format,  only  valid  for  scenarios  created  with  the  antenna  pattern,  shows  the  original  unnulled 
antenna  pattern,  the  nulled  antenna  pattern,  and  an  optimally  nulled  antenna  pattern.  This  plot 
is  repeated  for  each  weight  vector  in  the  results  file.  The  optimal  pattern  is  the  pattern  which 
would  have  been  obtained  if  the  weight  vector  had  been  solved  using  the  ideal  Cholesky  factor  and 
using  floating  point  Givens  rotations  to  solve  for  the  weight  vector.  The  second  plot  shows  the  S/N 
improvement  vs  time.  In  this  plot,  the  S/N  improvement  is  calculated  for  each  weight  vector  using 
the  appropriate  scenario  matrix,  and  then  the  S/N  is  plotted  vs  the  cycle  number. 

11.5.4  Simulation 

The  software  includes  a  bit-by-bit  simulation  of  the  hardware  to  allow  the  results  of  the  MUSE 
system  to  be  checked  for  accuracy.  Once  the  user  has  specified  the  addressing  information,  the 
snapshot  cycles  and  the  scenarios,  a  simulation  of  the  MUSE  may  be  performed  instead  of  running 
the  data  through  MUSE.  The  simulation  produces  the  same  type  of  results  file  as  discussed  above, 
and  therefore,  the  results  may  be  viewed  using  the  methods  discussed  above. 
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12.  SUMMARY 


This  report  has  described  the  complete  design  of  a  wafer  scale  adaptive  nulling  processor 
called  MUSE,  whose  individual  cells  carry  out  coordinate  rotations  using  the  CORDIC  algorithm. 
The  CORDIC  cells  are  almost  100  percent  utilized  and  communicate  with  one  another  in  a  simple 
systolic  fashion.  The  high  efficiency  and  modularity  of  the  algorithm  is  achieved  by  interleaving 
two  data  streams  travelling  in  opposite  directions,  and  depending  on  a  careful  choice  of  CORDIC 
latency  to  avoid  collisions. 

A  CORDIC  cell  including  memory  requires  approximately  54,000  CMOS  transistors  and  occupies 
a  5.5  mm  x  5.6  mm  rectangle.  It  has  been  clocked  at  12  MHz,  so  that  an  entire  system  of  96 
such  cells  can  update  a  64-element  weight  vector  for  300  observations  in  6.7  ms.  A  conventional 
computer  would  need  to  be  capable  of  2.88  Giga-ops  to  carry  out  the  same  task.  Simulations 
have  demonstrated  that  the  system  can  support  50  dB  of  S/N  improvement.  A  yield  analysis  has 
demonstrated  that  the  entire  MUSE  system  can  be  realized  on  a  single  wafer  if  cell  yield  averages 
70  percent.  Furthermore,  the  simple  systolic  communication  permits  using  two  or  more  wafers  if 
cell  yield  is  too  low. 

The  CORDIC  cells’  main  function  is  to  update  the  Cholesky  factor  of  the  correlation  matrix  of 
observed  interference  as  new  data  is  observed,  with  a  new  update  approximately  every  22  ns.  The 
solution  of  a  set  of  linear  equations  whose  coefficients  are  a  snapshot  of  said  Cholesky  factor  gives 
the  desired  nulling  weights.  The  CORDIC  cells  are  also  used,  in  an  interleaved  fashion,  to  solve  these 
linear  equations. 

A  short  chain  of  twrelve  CORDIC  cells,  packaged  as  chips  and  representing  one-eighth  of  the 
system  to  be  realized  on  a  wafer,  has  demonstrated  the  correct  operation  of  the  cells,  their  systolic 
communication,  system  timing,  and  S/N  improvement,  in  short,  all  aspects  of  MUSE  performance. 
The  wafer  floorplan  and  a  strategy  for  testing  the  wafer  as  it  is  restructured  are  completely  designed. 
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APPENDIX  A 

CREATION  OF  DATA  SETS  WITH  CONTROLLED  CONDITION  NUMBER 
AND  ACHIEVABLE  NULLING  PERFORMANCE 


Let  the  data  set  we  want  to  create  be  designated  as  U,  a  matrix  with  M  columns  and  N  rows. 
Its  columns  are  the  vectors  we  will  play  into  MUSE.  Any  nonsingular  matrix  may  be  expressed  in 
a  singular  value  decomposition 

U  =  EEV 

where  E  is  an  N  x  N  unitary  matrix,  V  is  an  M  x  M  unitary  matrix  and  E  is  an  N  x  M  rectangular 
diagonal  matrix  with  positive  elements,  E„  =  ax.  Also, 

<71  >  02  >  •  •  •  >  <7 ;v  >  0. 

In  what  follows,  the  matrix  V  will  have  only  a  very  small  role. 

The  correlation  matrix  for  this  data  set  is 


Ai 


R  =  EY?Eh  =  E 


Eh , 


Aw-i 


A,v  . 


with  A,-  =  of. 

One  of  the  quantities  we  want  to  control  is  the  condition  number  of  the  data  set.  The  condition 
number  is  We  can  always  scale  U  appropriately  so  that  ay  =  1.  Without  loss  of  generality, 
assume  that  we  have  done  so.  Therefore,  in  addition,  Aw  =  1. 

Once  a  weight  vector  W_  is  found,  by  whatever  means,  v  can  be  calculated,  which  is  the 
improvement  in  S/N  for  the  steering  vector  S.  Here  is  how  we  might  do  so  for  the  special  case  of  a 
sidelobe  canceller,  5  =  [0, 0, ... ,  1]*.  The  before  signal  power  is  ShS  =  1.  The  before  noise  power 
is  S_hRS.  The  after  signal  power  is  |SHU|2  =  |tuw|2-  The  after  noise  power  is  WhRW .  Thus  the 
S/N  improvement  for  a  general  weight  vector  is 


v  = 


ShRS 

WhRW 


|ww|2- 


The  optimum  weight  vector  is  the  solution  to 
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RW  =  kS 

where  k  is  any  constant.  Hence 

£hRS.  |2  , 

"opt  =  — rrl^l  'and 
wnk 

I'opt  =  ShRS 

In  this  last  equation,  w jy  depends  on  k.  We  can  choose  a  value  of  k  which  makes  come  out  to 
be  1.  In  terms  of  the  singular  value  decomposition  of  [/,  this  value  is 

t-  1 

Similarly 

N 

ShRS  =  ]T  \^Ni?crb 

i=  1 

Therefore 

V opt  =  (Y,  \eNi?o-2i)(Y  |ejv<|  2/crf ). 

1=1  1  =  1 

Now  it  is  possible  to  choose  U  such  that  its  singular  values  are  “controlled”  and  yet  consistent 
with  a  given  amount  of  S/N  improvement. 

Suppose  only  two  distinct  singular  values  are  used: 
c T{  —  a  for  1  <  i  <  Kand 
<jj  =  1  for  K  <  i  <  N. 

In  other  words,  we  are  modeling  K  equal  amplitude  and  orthogonal  jamming  sources.  The  we  must 
have 

"opt  =  {Y  |ew|V+  y  i^i2)(X>Nil2*-2+  Y  ie™i2)- 

i=l  i=K+ 1  i=l  i=K+i 

Abbreviate  \eNi\2  as  q  and  hence  £il/e+i  leJVt|2  as  (1  -  q).  Then 

"opt  =  (qv2  +  i  -  q){q<?~2  +  l-q). 
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Therefore 


9  ~9  + 


1  ~  vppt 


2  -  (cr2  +  a~2) 


=  0. 


The  approach  will  therefore  be  to  choose  the  desired  S/N  improvement  v^t  and  the  desired  condition 
number  a .  Then  we  can  solve  for  q,  which  tells  us  something  about  how  the  energy  in  the  last  row 
of  E  is  distributed.  E  can  then  be  computed  at  random  subject  only  to  this  constraint  and  to  the 
constraint  that  it  is  unitary. 

However,  not  every  possible  choice  of  v^t  and  cr  is  compatible,  q  must  be  real  and  positive. 
The  equation  for  q  will  give  complex  roots  unless 

1  ~  Vgpt  .  1 

2  -  (cr2  +  cr~2)  ~  A 


This  equation  gives  the  set  of  possible  scenarios  plotted  in  Figure  6-1.  The  boundary  = 
\  +  +  1  /g2)  is  also  plotted  in  Figure  6-1. 

To  construct  £*,  begin  by  constructing  its  last  row.  Without  loss  of  generality  this  last  row 
can  be  purely  real.  Begin  by  constructing  N  random  nonnegative  real  numbers  pj,  z  =  1 , 
Normalize  the  first  K  so  that  their  sum  is  q  and  normalize  the  last  1  —  K  so  that  their  sum  is 

(1-9)- 

K 

Tl  =  X>, 

1  =  1 

N 

T2  =  £  Pi 

i=K+ 1 

q 

Pi  *-  T?i 
1  -  q 

Pi  t2  Pi 

Now  set 


1  <  i  <  K 
K  +  l<i<N. 


em  <—  y/pi- 

Once  we  have  the  la<st  row  of  E ,  the  other  rows  can  be  determined  sequentially  by  a  Gram- 
Schmidt  orthogonalization  procedure.  For  each  new  row  we  generate  a  completely  random  complex 
Gaussian  precursor  row  Y 1  and  subtract  off  from  it  any  components  in  the  space  spanned  by  the 
rows  already  computed.  We  then  normalize  what  is  remaining  in  Y}  to  have  unit  length  (unless  its 
length  is  zero,  in  which  case  we  must  try  again  with  a  different  random  Y_ )  and  this  becomes  the 
next  row  of  E.  The  procedure  is  continued  until  we  have  all  N  rows  of  E. 
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When  we  have  E  and  £  (the  latter  being  a  diagonal  matrix  whose  diagonal  eleihents  are  cr 
and  1),  we  can  multiply  them  together  to  give  us  Uf  =  ££,  a  matrix  with  N  nontrivial  columns 
that  are  the  MUSE  inputs,  and  M  —  N  columns  which  are  all  zeroes.  We  may  interpret  the  first  K 
columns  of  Uf  as  “jammers”  turning  on  one  at  a  time,  each  with  energy  cr2  -  1,  and  the  columns 
K  -f  1  to  N  as  representing  noise  only. 

It  is  only  at  this  point  that  V ,  the  other  unitary  matrix,  is  needed.  V  allows  us  to  distribute 
the  noise  and  jammer  energy  defined  by  Uf  into  M  samples  in  any  way  we  like;  by  choosing  for 
V  we  can  choose  any  M  x  M  unitary  matrix,  without  affecting  the  two  constrained  properties  of 
U  =  f/'V\  its  condition  number,  and  the  available  S/N  improvement. 

One  interesting  possibility  is  to  choose  V  such  that  U  is  lower  triangular.  It  then  becomes  a 
data  set  that  is  its  own  Cholesky  matrix.  If  £T  is  played  into  MUSE,  we  expect  it  to  compute  L 
which  is  approximately  equal  to  that  Cholesky  matrix. 
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APPENDIX  B 
TESTING  A  MUSE  SUBSET 


The  choice  of  N  =  64  substantially  impacted  the  design  of  the  CORDIC  cell  for  MUSE.  One 
impact  is  due  to  the  relationship  between  N  and  latency  for  latency  controlled  interleaving.  A 
second  impact  is  in  the  size  of  memories  built  into  the  chip.  For  these  reasons,  the  design  of  the 
cell  dictates  the  system  timing. 

For  various  reasons,  it  may  be  desirable  to  use  the  CORDIC  cell  designed  for  N  =  64  to  realize 
an  adaptive  nulling  system  for  smaller  N.  In  the  body  of  this  report  two  such  requirements  are 
mentioned:  first,  the  desire  to  verify  the  operation  of  the  MITE  using  CORDIC  chips  that  were 
fabricated  to  prove  out  the  cell  design;  second,  the  testing  of  the  wafer  during  restructuring,  as 
additional  supercells  are  connected  into  the  working  system.  For  any  even  N9  less  than  N  the  MUSE 
design  can  be  used  for  adaptive  nulling  using  3iV' /2  CORDIC  cells  of  the  type  already  optimized  for 
N  =  64.  The  use  will  not  be  efficient. 

The  approach  we  take  is  to  “imbed”  the  length  Nf  nulling  problem  into  the  larger  length  N 
problem.  Figure  B-l  illustrates  how  an  Nf  —  4  input  vector  is  imbedded  in  the  middle  of  an  N  =  10 
input  vector,  and  how  an  N9  x  N '  Choleskv  matrix  is  therefore  imbedded  in  the  N  x  N  Cholesky 
matrix  developed  from  it.  The  three  leading  elements  of  the  input  vector  are  forced  to  zero  and 
the  last  three  elements  may  be  set  to  any  value — they  are  “don't  care”  elements  indicated  by  d’s  . 
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Figure  B-I.  4x4  adaptive  problem  imbedded  in  10  x  10  framework . 

If  the  ten-element  vectors  of  Figure  B-l  are  fed  into  a  mini-MUSE  made  of  five  supercells,  the 
first  three  supercells  in  the  five  supercell  chain  (supercells  0,  1,  and  2)  will  have  a  negligible  effect 
on  the  input  vectors  as  they  are  passed  forward  since  their  function  is  to  zero  out  the  elements 
which  are  already  zero.  Therefore  it  is  easy  to  predict  the  form  of  the  tacked-on  vector  that  will 
be  the  input  to  supercell  3,  namely  [xi,X2,X3,X4 ,  dTd,  d]. 

Therefore  we  can  provide  only  the  last  two  supercells  (in  general  the  last  N9/ 2  supercells),  and 
feed  in  the  predicted  vector  preceded  with  the  proper  number  of  blank  microcycles .  Assuming  that 
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the  correct  D,  Q,  and  ZERO  controls  are  provided,  the  two  supercell  MUSE  will  operate  exactly  as 
if  it  were  the  last  two  supercells  of  the  five  supercell  MUSE,  and  the  imbedded  Cholesky  factor  will 
be  computed.  Signals  can  also  be  provided  which  initiate  the  Q-operation  and  the  Q-op  controls 
will  be  emitted  from  supercell  3,  in  its  backward  direction,  so  that  the  weights  can  be  computed. 

The  example  with  two  actual  supercells  from  a  chain  of  five  supercells  easily  and  accurately 
extends  to  the  case  of  any  smaller  subset  of  a  chain  of  32  supercells.  There  are  three  points  to 
remember: 

•  The  actual  supercells  must  be  consecutive  and  must  include  supercell  31. 

•  The  timing  must  use  67  microcycles  per  macrocycle  and  the  input  vector, 
direction ,  Q,  and  ZERO  must  appear  as  they  would  appear  if  the  nonexistent 
supercells  were  present. 

•  The  leading  supercell  actually  present  must  be  provided  with  some  number  of 
slave  Q-operations  to  prompt  it  to  initiate  a  master  Q-operation  at  the  correct 
time. 


/ 
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